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I. Program Objectives 


In this program we proposed to perform a series of spectroscopic studies, including data analysis and 
modeling, of star formation regions using an ensemble of archival space-based data from the Infrared 
Space Observatory’s Long Wavelength Spectrometer and Short Wavelength Spectrometer, and to 
take advantage of other spectroscopic databases including the first results from SIRTF. Our empha- 
sis has been on star formation in external, bright IR galaxies, but other areas of research have in- 
cluded young, low or high mass pre-main sequence stars in star formation regions, and the galactic 
center. The OH lines in the far infrared were proposed as one key focus of this inquiry because the 
Principal Investigator (H. Smith) had a full set of OH IR lines from ISO observations. It was 
planned that during the proposed 2-1/2 year timeframe of the proposal other data (including perhaps 
from SIRTF) would become available, and we intended to be responsive to these and other such spec- 
troscopic data sets. 

The program has the following goals: 

1) Refine the data analysis of ISO observations to obtain deeper and better SNR results on selected 
sources. The ISO data itself underwent pipeline 10 reductions in early 2001, and the more 
“hands-on data reduction packages” have been released. The ISO Fabry-Perot database is par- 
ticularly sensitive to noise and can have slight calibration errors; improvements are anticipated. 
We plan to build on these deep analysis tools and contribute to their development. 

2) Model the atomic and molecular line shapes, in particular the OH lines, using revised montecarlo 
techniques developed by the Submillimeter Wave Astronomy Satellite (SWAS) team at the Cen- 
ter for Astrophysics. 

3) Use newly acquired space-based SIRTF or SOFIA spectroscopic data as they become available, 
and contribute to these observing programs as appropriate. 

4) Attend scientific meetings and workshops. 

5) Conduct E&PO activities, especially as related to infrared astrophysics and/or spectroscopy. 

The program since its inception has been extremely successful: We presented or submitted more 
than about a dozen articles, sponsored a workshop, and participated in scientific conferences. The 
molecular line montecarlo code has undergone important revisions for infrared lines and has made 
significant predictions. First Spitzer data are available and we have made useful contributions to 
pending SIRTF science observations. E&PO programs have been active, including the exhibit at 
Boston’s Museum of Science. We anticipate progress to continue in all areas. Detailed descriptions 
of all these results follow below. 


II. Progress Report 

We made excellent progress during the past year; however we note that funding for the year’s activi- 
ties did not arrive until last week — nearly 11 months late . As a result the progress was considerably 
below our hopes. In the next year we plan to complete the tasks unfinished due to this funding de- 
lay. 

A. ISO Data Analysis and Modeling 

Last year we submitted a paper entitled “The Far Infrared Lines of OH as Molecular Cloud Diagnos- 
tics” to the 2 nd Maryland Conference on Far IR Astronomy, 2002, and also presented an improved 
version of “The Infrared Lines of OH: Diagnostics of Molecular Cloud Compositions in Infrared 
Bright Galaxies,” to the Proceedings of the Workshop of the Neutral ISM in Starburst Galaxies, Mar- 
strand, Sweden, 2003. This analysis was based on the SWAS montecarlo code modeling of the OH 
lines in galaxies observed by ISO (Figure 1). Since that meeting last spring considerable effort has 
been put into improving the montecarlo code. A group of European astronomers, including Dr. Pepe 
Chernicharo and Prof. Eduardo Gonzdlez-Alfonso, had been performing montecarlo modeling of other 
molecules seen in ISO galaxies We used portions of this grant to bring Prof. Gonzalez to Cambridge 
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for an intensive working visit. The result of that work is the paper in Appendix C on Arp220, now 
submitted to ApJ. 

In a second major paper on the IR spectroscopy of galaxies, we are about to submit a paper on the 
spectroscopy on NGC 1068, the archetypal Seyfert 2 galaxy (Figure 2). One of the dramatic results 
in our work is the discovery that ALL of the OH lines seen there are seen in emission, making this 
source unique. A preprint is attached in Appendix C. 

B. Other Science Results 

SIRTF : SIRTF science development was the other major component of this year’s research. This 
program supported the development of five “Early Release Objects” for SIRTF observations on which 
Dr. Smith was Principal Investigator or Co-Investigator, and another 5 proposals for GO time. The 
early release program is designed to rapidly present to the public and the scientific community some 
exciting results from Spitzer in the first months of its operation. The Spitzer instrument and science 
teams submitted proposals for ERO objects, and a competitive selection process narrowed these 
down to a small group with exciting science and realistic observational parameters. This grant sup- 
ported Dr. Smith’s participation in the ERO process, including developing science goals, identifying 
key objects for observation, and developing the detailed AOR (observing formulae) to be used by the 
instruments for mapping, integrating, etc. 

Appendix B lists the Spitzer programs involved. 

C. Advising Students and Postdocs 

During this year we had one graduate student, Idan Ginsburg, who participated in the modeling of 
the continuum of galaxies using the montecarlo code and who wrote a C++ script to facilitate this 
work. We also supervised the visit of Dr. Eduardo Gonzdlez-Alfonso, who likewise worked on the gal- 
axy modeling (reprints attached). 

D. E&PO Activities 

Diming this period the PI gave several talks to elementary school groups, as well as to a summer 
camp. 


III. Program Plans 

During the final year of this program, much of which has been delayed because funding arrived 

eleven months late, we intend to bring to closure the projects discussed to date: 

1) Reduction of the ISO galaxy data set using the Chernicharo/Gonzdlez algorithms; 

2) Modeling the OH lines with the corrected and improved montecarlo code; 

3) Tightening the arguments made in our earlier article -- that there are three types of active/IRB 
galaxies as typified by the kinds of molecular clouds that dominate their FIR emission: large dif- 
fuse clouds as in normal IRBs; smaller, more dense clouds as in Arp 220 and other mixed AGN; 
and very dense small cloudlets as seen in the Seyfert galaxy NGC 1068; 

4) Spitzer activities: We will continue to be active in E&PO and GO, including an aggressive 
IRAC/Spitzer E&PO effort. 


IV. Reportable Patents/New Technology 


There have been no reportable patents/new technology made under this grant. 



-18 


12 


1/2— 3/2 


1 

I T7 TT 

4 22“ 3 13 

[0 1], 

H O 


1 n i/2 n 3/2 

3/2-3/2 


n 2 u 

1 33.-2 20 

oh n 3/2 

9/2— 7/2 

oh n 1/2 

7/2— 5/2 



(4,3).-(3,3). 


(4.2).-(3.2). 


(4.1).-(3.1). 


(4,0)— (3,0) 

oh n 3/z 

NH S ✓ 

5/2— 3/2 

(4.3),— (3,3), 

1 

(4,2),— (3,2), 

, 1 

(4.1),— (3.1), 



Figure 1 









\(fi m) 


Fig. 1. — The ISO-LWS spectrum of NGC1068. The long-dashed line, which is fitting the 
data (reduced x 2 =0.40, corresponding to a probability of P>98%) represents the sum of the 
three different grey-body functions: dotted line: main component at T=32K; dotted-dashed 
line: cold component at T=20K; short-dashed line: warm component at T=65K. 


Figure 2 
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1) IRAC Images of an Infrared Bright Galaxy: M82 

2) IRAC/IRS Early Release Observations of RCW 108 

3) Observations of Interacting Galaxies 

4) SIRTF Observations of the Galactic Center 

5) Extremely Red Background Objects in the IRAC Survey 

Spitzer GO Proposals 

6) The Spitzer Interacting Galaxy Sample 

7) The Low Redshift 3CRR Sources 

8) Probing the Dynamical Causes of the Morphological Butcher-Oemler Effect 

9) IRS Spectroscopy of a Complete Sample of Active and Starburst Galaxies 

10) Neon, Sulfur and Discrepant Abundance Gradients in M33 
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1) The Infrared Lines of OH: Diagnostics of Molecular Cloud Conditions in Infrared Bright 
Galaxies. 

2) The Far-Infrared Spectrum of Arp 220 

3) The Far-Infrared Emission Line and Continuum Spectrum of the Seyfert Galaxy NGC 
1068 


**TITLE** 

ASP Conference Senes, Vol. “VOLUME**, **PUBLICATION YEAR** 
“EDITORS** 


The Infrared Lines of OH: Diagnostics of Molecular Cloud 
Conditions in Infrared Bright Galaxies 


Howard A. Smith(l), Eduardo Gonzalez- Alfonso(l, 2), Jackie Fischer(3), 
Matthew L.N. Ashby (1), Chris Dudley(3), and Luigi Spinoglio(4) 

(l)Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 
Cambridge, MA 02138, (2)Universidad de Alcala, Spain. (3)Naval 
Research Laboratory, Washington, DC 20375, (4)Istituto di Fisica dello 
Spazio Interplanetario, CNR, Italy 

Abstract. ISO detected 9 of the 15 far-IR lines of OH that arise be- 
tween the molecule’s lowest 8 rotational levels. The lines were seen, some- 
times in emission and sometimes in absorption, in a wide variety of galac- 
tic and extragalactic objects ranging from AGB stars to molecular clouds 
in AGN and LIRGs. This set of OH lines provides an important diagnos- 
tic for many reasons: they span a wide wavelength range (from 28.9pm to 
163.2pm); the transitions have fast radiative rates; the abundance of OH 
is relatively high; the IR continuum plays an important role as a pump; 
the contribution from shocks is relatively minor; and, not least, the strong 
cm-wave radiation from OH allows comparison with VLBI datasets. OH 
plays a key role in the chemistry of the ISM. We have used montecarlo 
and non-local radiative transfer codes to analyze the spectra of a number 
of LIRGs with excellent success., Other detected lines in the ISO spec- 
tra, including those from H2O, CH, NH3, CO, [CII] and [01], are also 
modeled by the codes, as is the continuum. 


1. Characteristics of clouds in LIRGs as inferred from their OH lines 

ISO-LWS detected OH in emission and/or absorption in: Arp 220, Cen A, 
IRAS17208-0014, IRAS20100-4156, M82, Mkn 231, Mkn 273, NGC 253, NGC 
891, NGC 1068, NGC 1614, NGC 3690A, NGC 4945, NGC 7469, and 3Zw35. 
Perhaps the most striking result from the observed set of extragalactic lines is 
their wide range of behavior, even though all arise in LIRGs with either active 
star formation, an active nucleus, or perhaps both. In Arp 220, for example, 
every OH line is seen in absorption except the longest wavelength, 163pm line 
which is seen weakly in emission (Gonzalez- Alfonso et al., 2004; Fischer et al., 
’’Proc.ISO to the Peaks Workshop,” 1998). By contrast, in NGC 1068 every 
detected line is seen in emission, even the strong 119pm transition between 
ground and the first excited state (Spinoglio et al., 2004). NGC 253, a nearby 
starburst, has some OH lines in emission and some in absorption (Bradford et 
al., ’’Proc.The Universe as Seen by ISO,” 1999; Smith et al., ”Proc. 2nd Md. 
Conference on FIR Astronomy,” 2004), while in M82, the infrared bright, pro- 
totype starburst galaxy, the lines’ equivalent widths are so small that even with 
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Figure 1. The full ISO-LWS spectra of NGC 253 (strong [CII] ) and 
Arp220 (strong water), with the prominent features marked. The con- 
tinua and all the features are well modeled (from Gonzalez- Alfonso et 
al., 2004). 


our high signal-to- noise ratio only the 119^/m line has been conclusively seen at 
all, in absorption (Colbert et al., Ap.J., 511, 721, 1999). These differences in 
the OH lines reflect important differences in the molecular clouds that dominate 
in each of these various types of galaxies, and OH is a useful tool to distinguish 
these cloud properties - and the galaxies’ morphologies. 

We used two codes, a MonteCarlo code developed for SWAS (Ashby et 
al., Ap.J.Lett, 539, L119, 2000) and a non-local code (Gonzalez-Alfonso and 
Chernicharo, A&A, 322, 938,1997), as well as DUSTY (for the continuum), to 
calculate and cross-check line intensities and profiles in spherical clouds. The 
source is divided into a series of independent, concentric shells with adjustable 
temperatures, densities, abundances, etc.; microturbulence and/or radial veloc- 
ity gradients are also included. The results confirm that the OH lines are often 
very optically thick (exceeding 100 in many cases). Since the critical densi- 
ties for OH are high, above about 10 8 cm -3 , the level populations are primarily 
determined by radiative effects. 

From the analysis of this diverse set of sources and line strengths, several 
useful preliminary generalizations may be drawn. We tentatively conclude that 
galaxies can be grouped into three broad categories based solely on their OH 
line behavior. AGN: For NGC 1068, we model the clouds as being small, only 
a few pc in size, and dense. ULIGS: In Arp 220 the 34/xm absorption confirms 
the OH maser radiative pumping mechanism (Skinner et al., Nature, 370, 117, 
1997), and we find high OH column densities of about 2 x 10 17 cm~ 2 . LIRGs: For 
galaxies like NGC 253 and Mkn 231, which have differing combinations of both 
emission and absorption lines, we typically find that the clouds are much larger, 
10-100pc, and less dense. We have made substantial progress in the analyses, 
and future velocity-resolved spectra will enable even more refinements. Sensitive 
OH measurements should help categorize those sources whose morphology is 
otherwise ambiguous. 
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ABSTRACT 

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 shows absorption 
in molecular lines of OH, H2O, CH, NH, and NH3, as well as in the [0 I] 63 line and emission 
in the [C II] 158 pm line. We have modeled the continuum and the emission/absorption of all 
observed features by means of a non-local radiative transfer code. The continuum from 25 to 1300 
pm is modeled AS A WARM (106 K) NUCLEAR REGION THAT IS OPTICALLY 
THICK IN THE FAR-INFRARED, attenuated by an extended region (size 2") that is 
heated mainly through absorption of nuclear infrared radiation. The molecular absorption in the 
nuclear region is characterized by high excitation due to the high infrared radiation density. The 
OH column densities are high toward the nucleus (2 - 6 x 10 17 cm -2 ) and the extended region 

2 x 10 17 cm'' 2 ). The H2O column density is also high toward the nucleus (2 — 10 x 10 17 
cm* 2 ) and lower in the extended region. The column densities in a halo that accounts for the 
absorption by the lowest lying levels are similar to what are found in the diffuse clouds toward 
the star forming regions in the Sgr B2 molecular cloud complex near the Galactic 
Center. Most notable are the high column densities found for NH and NH3 toward the nucleus, 
with values of — 1.5 x 10 16 cm" 2 and ~ 3 x 10 16 cm" 2 , respectively, whereas the NH2 column 
density is lower than ~ 2 x 10 15 cm -2 . A combination of PDRs in the extended region and 
hot cores with enhanced H 2 O photodissociation and a possible shock contribution 
in the nuclei may explain the relative column densities of OH and H2O, whereas the nitrogen 
chemistry may be strongly affected by cosmic ray ionization. The [C II] 158 /im line is well 
reproduced by our models and its “ deficit ” relative to the CII/FIR ratio in normal and starburst 
galaxies is suggested to be mainly a consequence of the dominant non-PDR component of far- 
infrared radiation, ALTHOUGH OUR MODELS ALONE CANNOT RULE OUT EXTINCTION 
EFFECTS IN THE NUCLEI . 
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Subject headings : galaxies: abundances — galaxies: individual (Arp 220) — galaxies: ISM — galaxies: 
starburst — infrared: galaxies — radiative transfer 



1. Introduction 

With a redshift of z = 0.018, Arp 220 (IC 
4553 / 4 ) is the nearest and one of the best 
studied ultraluminous infrared galaxies (ULIRGs). 
The tails observed in the optical, together with the 
double highly-obscured and compact (0"3) nuclei 
observed in the near and mid-infrared, as well as 
in the millimeter, strongly suggest that the enor- 
mous luminosity of Arp 220, ~ 10 12 L©, is the 
result of galactic merging. Nevertheless, the con- 
crete physical process responsible is still a matter 
of debate: the proposed sources are hidden active 
nuclei and/or bursts of star formation. 

Molecular observations of Arp 220 provide 
unique clues to the physical and chemical pro- 
cesses ocurring in the nuclei and their surround- 
ings. In the millimeter region, CO observations 
have been carried out with increasingly high an- 
gular resolution (Radford, Solomon, & Downes 
1991; Scoville et al. 1991; Okumura et al. 1994; 
Scoville, Yun, & Bryant 1997; Sakamoto et al. 
1999). In particular they have shown that, on the 
one hand, CO emission arises from a region signif- 
icantly more extended than the nuclei (~ 3 - 4"), 
and on the other hand that the CO(2-l) to CO(l- 
0) intensity ratio is lower than 1, thus suggest- 
ing that CO mainly traces low density regions 
(< 10 1 2 3 cm" 3 ). Observations of molecules with 
high dipole moment, like CS and HCN, have re- 
vealed that the fraction of molecular gas contained 
in dense clouds (n(H 2 )> 10 4 cm -3 ) is much larger 
than in normal galaxies, yielding ~ 10 10 M© of 
dense gas (Solomon, Radford, & Downes 1990; 
Solomon, Downes, & Radford 1992). Radford et 
al. (1991b) found that HCN(l-O) and HCO+ (1-0) 
peak strongly toward the nuclei, but also show 
low-level extended emission. More recently, Aalto 
et al. (2002) have detected emission from the high 
density tracers HNC and CN, and the relatively 
low HCN/CN and HCN/HNC intensity ratios 
were attributed to widespread PDR chemistry. 


1 Based on observations with the Infrared Space Obser- 
vatory, an ESA project with instruments funded by ESA 
Member States (especially the principal investigator coun- 
tries: France, Germany, Netherlands, and the United King- 
dom) and with the participation of ISAS and NASA. 

2 Visiting Astronomer, Harvard- Smithsonian Center for 
Astrophysics, 60 Garden Street, Cambridge, MA 02138. 

3 CSIC, IEM, Dpto. Astrofisica Molecular e Infrarroja, 

Serrano 123, E-28006 Madrid, Spain. 


The launch of the Infrared Space Observatory 
(ISO) opened a new window for the study of the 
physical and chemical properties of ultraluminous 
infrared galaxies. Despite the lack of angular and 
spectral resolution, the observations of Arp 220’s 
far-infrared spectrum from 40 to 200 /xm (Fischer 
et al. 1997, 1999) and of a number of individ- 
ual lines in the SWS range (Sturm et al. 1996) 
provided new insights in our understanding of the 
ionic, atomic and molecular content of the galaxy. 
These wavelength regions are of great interest, be- 
cause the bulk of the enormous luminosity is emit- 
ted in the far-infrared, and also because they con- 
tain lines of interesting molecular species. Skinner 
et al. (1997) reported the detection of the 35 /zm 
OH line in Arp 220. Fischer et al. (1997, 1999) 
found that the far-infrared molecular absorp- 
tion lines of OH, H 2 0, CH, and NH 3 are 
significantly stronger in Arp 220 than in 
less luminous infrared-bright galaxies while 
the fine structure lines from ionic species 
are, to the contrary, extremely weak. Luh- 
man et al. (1998, 2003) found that, rela- 
tive to the far-infrared luminosity, the [C 
II] in ULIRGs is typically nearly an order of 
magnitude weaker than in lower luminosity 
infrared-bright galaxies. Sturm et al. (1996) 
reported the detection of two ortho-H 2 pure rota- 
tional lines, indicating that high masses of gas are 
subject to PDR conditions and/or shock activity. 

The physical and chemical processes that ac- 
count for the rich molecular far-infrared spec- 
trum of Arp 220 can be better understood if quan- 
titative values of the column densities of the above 
species, as well as their excitation conditions, are 
estimated. The presence of OH and, to some ex- 
tent, of H 2 0, may be indicative of PDR and/or 
diffuse interstellar cloud chemistry, and their col- 
umn densities potentially give an estimate of the 
UV field in the source. On the other hand, large 
amounts of H 2 0 are produced in non-dissociative 
shocks (e.g. Cernicharo et al. 1999), where the 
OH abundance is also enhanced (Watson et al. 
1985). The OH abundance is expected to be 
generally higher than that of H 2 0 in fast disso- 
ciative shocks (Neufeld & Dalgarno 1989). H 2 0 
ice in grain mantles may also efficiently return 
to the gas phase through sublimation of mantles 
in “hot core” regions. But whatever the char- 
acteristics of the regions producing the observed 
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molecular features, the lines under study lie at 
wavelengths where the enormous infrared contin- 
uum flux approaches its maximum, and the molec- 
ular excitation of high-dipole species should be 
strongly affected by absorption of continuum radi- 
ation. Hence any reliable estimation of molecular 
column densities require accurate models for the 
dust emission. Unfortunately ISO’s lack of angu- 
lar resolution forces us to rely on plausibility argu- 
ments in our assumptions about the regions where 
the different lines are formed; some of these are 
based on general requirements of excitation, and 
others on conclusions from observations of galac- 
tic sources. The main goal of this work is thus 
to shed light on the physical and chemical pro- 
cesses in Arp 220, based on detailed model fits of 
its continuum and far-infrared molecular/atomic 
line absorption and emission spectrum. We adopt 
a distance to Arp 220 of 72 Mpc (projected linear 
scale of 350 pc/arcsec, Graham et al. 1990). In 
section 2 we present the ISO/LWS observations; 
in section 3 we discuss the line identifications; 
section 4 is a discussion of the models for the 
continuum emission; section 5 presents the models 
for the molecular and atomic species; in section 6 
we discuss the implications of the radiative trans- 
fer models, and section 7 summarizes our main 
results. 

2. Observations 

The full 43-197 pm spectrum of Arp 220, 
obtained with the LWS spectrometer (Clegg et 
al. 1996) on board ISO (Kessler et al. 1996) 
(TDT— 27800202), was presented by Fischer et 
al. (1997, 1999). The grating spectral resolution 
is £s 0.3 pm in the 43-93 pm interval (detectors 
SW1-SW5), and of ?^0.6 in the 80-197 pm in- 
terval (detectors LW1-LW5), corresponding to 
Av > 10 3 kms" 1 . The lines are thus not resolved 
in velocity. The beam size of « 80" ensures that 
all the continuum and line emission/absorption 
from Arp 220 (CO size < 4", Scoville et al. 1997, 
hereafter SYB97) lie within the ISO aperture. 

The data we present here was reduced using 
version 9.1 of the Off Line Processing (OLP) 
Pipeline system which we found to produce 
a higher signal-to-noise spectrum than OLP 
10.1. However, we adopted the continuum 
correction given by OLP version 10.1, which 


typically gives absolute responsivity correc- 
tions with uncertainty factors ~ 3 times 
lower than are produced by version 9 (Tim 
Grundy, private communication) . In order to 
obtain a smooth spectrum throughout the whole 
LWS range, shifting scale factors lower than 10% 
were applied to several detectors, except for de- 
tector SW1 (43-64 pm) for which the correction 
factor was of 15%. Thus we attribute a conserva- 
tive uncertainty of 15% to the overall continuum 
level. 

The LWS spectrum of Arp 220 is pre- 
sented in Figure 1 together with identifica- 
tions of the most prominent lines. Owing 
to transient effects, the fluxes of weak lines as ob- 
served in the forward and reverse scans were found 
to differ significantly in some wavelength ranges. 
In these cases, the scan that showed higher signal- 
to-noise /bf ratio, or the average if the ratio was 
similar, was used. Nevertheless, the reverse scan 
was selected if the line appeared close to the upper 
end of the detector, because transient effects are 
expected to be minimal there. In wavelength re- 
gions where two detectors’ responses overlap, line 
fluxes were generally found to be consistent. The 
only exception was the H 2 0 3 22 - 2 U line at 90 
pm, which showed in LW1 a flux 60% weaker than 
in SW5. We adopted here the SW5 spectrum, but 
the flux of the above H 2 0 line should be consid- 
ered highly uncertain. The subtraction of a base- 
line (see Fig. 1) added additional uncertainty to 
the line fluxes, particularly in cases of broad fea- 
tures presumably composed of several lines. With 
the exception of the H 2 0 3 22 - 2n line, we esti- 
mate a line flux uncertainty generally lower than 
35%. 

3. General results 

OH and H 2 O: The FIR spectrum of Arp 
220 is dominated by lines of OH and H 2 0 
that are observed in absorption, with the ex- 
ception of the OH n 1/2 3/2 - 1/2 emission line at 
163.3 pm. Figure 2 shows the level diagram of OH, 
ortho-H 2 0 and para-H 2 0, and indicates the lines 
detected in Arp 220. Lines with very different ex- 
citation requirements are observed throughout the 
spectrum. The OH lines n 3 / 2 J = 9/2 - 7/2 at 65 
pm and U X / 2 J = 7/2 - 5/2 at 71 pm have lower 
levels at 290 K and 415 K above the ground state, 
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Fig. 1. — ISO/LWS spectrum of Arp 220, where the most prominent line features are identified (see text) 
The grey line shows the adopted baseline (continuum level). 
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respectively, whereas the H 2 O lines 4 32 - 3 2 i (59 
/zm) and 4 22 - 3 i3 (58 /zm) have lower levels at 
305 K and 205 K. Strong absorption is also ob- 
served in the OH ground state lines at 53, 79, and 
119 /zm, as well as in the H 2 0 lowest-lying line 
at 179 /zm. This wide range of excitation sug- 
gests that several regions with different physical 
conditions are contributing to the observed fea- 
tures (Fischer et al. 1999), and one of the goals of 
this work is to provide a reasonable estimate of the 
nature of these regions and their relative contribu- 
tions to the spectrum. On the other hand, several 
of the lines have complex shapes, with evidence 
of shoulders suggestive of weaker secondary lines. 
In particular, the ground state OH 119 /zm line 
shows a redshifted “shoulder”, which is detected 
in both the forward and reverse scans, although 
with somewhat different strengths. It could be 
attributed to the 18 OH ri 3/2 5/2 — 3/2 line, al- 
though contamination by other species such as 
CH + cannot be ruled out. Also, the redshifted 
“wing” of the H 2 0 2 12 - loi line at 179 /zm, at- 
tributed in Fig. 1 to the H 2 0 2 2 i - 2 12 line, could 
also be contaminated by H^O, CH (see Fig. 7), 
andH 3 0+ 

CH and NH: The spectrum contains lines from 
other molecular species: like CH at 149 /zm, NH 3 
at 125, 166 and 170 /zm and, very interestingly, 
strong absorptions at 102 and 153.2 /zm that have 
been identified as NH in Fig 1. Evidence for 
the latter identifications is strengthened be- 
cause of the presence of weak line-like features at 
155.74 /zm, and marginally at 151.53 /zm, which 
would correspond to the NH 2i — li and 2\ — lo 
lines, respectively (see also Fig. 3). Conceiv- 
ably, the line absorptions at 102 and 153.2 /zm 
could be severely contaminated by other species, 
like C 3 , H^ 8 0, NH 3 , and even OH 4- . C 3 has a 
strong transition at 153.3 /zm, but its contribu- 
tion is expected to be minimal due to the lack of 
detection of other adjacent C 3 lines (in particular 
at 154.86 /zm). The absorption at 102 /zm may 
be contaminated by the H^O 2 20 — I 11 line (just 
at 102.0 /zm), but since the H^O should not be 
as strong as the corresponding adjacent line of the 
main isotope H^O, we regard this identification 
as also unlikely. Some contribution of NH 3 lines 
to the 102 /zm feature might be expected, but they 
are somewhat shifted in wavelength (they lie be- 
tween 101.5 and 101.7 /zm). Both features could be 


contaminated to some extent by OH + , with strong 
lines at 101.70 and 101.92 /zm ( Nj = 3 3 — 2 2 and 
34 - 2 3 ) and at 153.0 /zm (2 3 - 1 2 ), but the strong 
absorption at 153.2 /zm absorption could never be 
explained by OH + alone. Therefore, despite the 
possible contamination from other molecules, NH 
is probably responsible for most, if not all, of 
the observed 153.2 /zm absorption (see also Fig. 
7). Thus, although the definitive assignment 
to NH should await confirmation with higher 
spectral resolution observations of the lowest-lying 
NH transitions at ~ 10 3 GHz, we conclude that 
ISO observations strongly support its detection in 
Arp 220, and advance a model of the observed 
absorption (section 5) that can be useful to di- 
rect future observational and theoretical studies. 
If confirmed, this detection is the first extragalac- 
tic detection of NH, which has been previously de- 
tected only in galactic diffuse clouds through elec- 
tronic transitions (Meyer &: Roth 1991; Crawford 
& Williams 1997) and, interestingly, toward Sgr 
B2 in the Galactic Center via the same transitions 
detected in Arp 220 (Cernicharo, Goicoechea, & 
Caux 2000; Goicoechea, Rodriguez-Fernandez, & 
Cernicharo 2004, hereafter GRC04). 

NH$ : The spectrum of Arp 220 around 125 and 
170 /zm, shown in Figure 4, strongly supports the 
identification of NH 3 . The shape of the 165.7 /zm 
feature indicates transient effects, but the line is 
detected in both the forward and reverse scans. 
Some H 2 0 lines may contribute to the observed 
125 and 127 /zm features, but they are shifted in 
wavelength relative to the strongest absorption. 
This is the first extragalactic detection of infrared 
NH 3 lines. With the detection of NH and 
NH 3 we might expect to detect NH 2 , but 
there is no evidence for its strongest expected 
lines at 159.5, 117.8 and 104.9 /zm. 

CO?: An apparent emission line, detected 
in both the forward and reverse scans, is 
present at 173.7 /zm in Fig. 1 . It coincides 
rather well with the expected position of the CO 
J = 15 — 14 line at 173.63 /zm. This identi- 
fication cannot be confirmed by the detec- 
tion of other expected CO lines because the 
CO J — 14-13 line at 186.0 /zm lies at the 
noisy edge of the LW5 detector and the CO 
J — 16—15 line at 162.8 /zm is blended with 
the OH 163.3 /zm line. The higher J lines are 
expected to be too weak to be detectable, given 
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OH 


ortho -H 2 0 


para-H 2 0 



Fig. 2. — Energy level diagrams of OH and H2O (ortho and para). Rotational levels with energies above 
the ground state up to 620 K are shown; the lines detected in Arp 220 are indicated with arrows and their 
wavelengths are in /xm. OH A-doubling is ignored because the A-doublets are not resolved with the ISO 
grating resolution. 




Fig. 4. — NH3 lines around 125, 127, 166, and 170 /im. The grey line shows the adopted continuum level. 
The 125 and 127 /xm features could be partially contaminated by the labelled H 2 0 lines. 
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the observed strength of the 173.7 /zm feature. 

[0 1] and [C II]: The ISO spectrum also shows 
the [O I] 63.2 /zm line in absorption and the [C 
II] 157.7 /im line in emission. The [0 I] 145.5 /im 
line is not detected. These lines, as observed in 
Arp 220 and other ULIRGs, have been discussed 
elsewhere (Fischer et al. 1997, 1999; Luhman et 
al. 1998, 2003). In section 5 we present a simple 
model of the Arp 220 spectrum that may shed 
some light on the peculiar behavior of these lines 
in Arp 220. 

In summary, the far-infrared spectrum of Arp 
220 shows molecular lines of OH, H 2 O, CH, NH, 
and NH 3 . The atomic lines of [0 I] at 63 /zm and 
[C II] at 158 /zm are also detected. Lines of other 
species, like CO, HsO -1 ", CH + , and OH + , could 
also contaminate the observed features, but our 
limited spectral resolution prevents the possibility 
of unambiguous detection. Only the [C II] 158 /zm 
and the OH Ui/ 2 3/2- 1/2 163 /zm lines are clearly 
observed in emission. Lines from ions that would 
trace H II regions and/or an AGN are absent. 

4. Models for the continuum 

Figure 1 shows that the continuum peaks 
around ~ 40 - 50 /im. The bulk of the con- 
tinuum from Arp 220 is emitted by heated dust 
grains. At 1.3 millimeter wavelengths, Sakamoto 
et al. (1999) showed that the continuum arises 
almost exclusively from the nuclei, with an equiv- 
alent size of ~ 0!'4. On the other hand, Soifer et 
al. (1999, hereafter S99) have shown that the two 
nuclei also account for essentially all the contin- 
uum at 25 /im. Combining both observations, S99 
proposed two alternative scenarios to explain the 
continuum emission of Arp 220 from far-infrared 
to millimeter wavelengths. Our models of the 
continuum emission are entirely based on these 
scenarios, which we have examined and refined 
quantitatively on the basis of our ISO 45-200 /im 
spectrum. 

Model Si : In the first scenario (hereafter Si), it 
is assumed that the emission from the nuclei is not 
significantly attenuated at 25 fim by foreground 
material. We have simulated the emission from 
the nuclei as arising from a single nucleus with 
effective size of 0!'41. With an effective dust tem- 
perature of 85 K, and optically thick emission in 
the submillimeter, the requirement that the fluxes 


at 25 and 1300 /im arise from the nuclei is ful- 
filled. However, the emission at 60-100 /zm, as 
well as the total luminosity from the galaxy, are 
then underestimated and a more extended region 
(hereafter ER) must be invoked to account for the 
remaining flux. We identify this surrounding en- 
vironment with the extended emission observed in 
CO, HCN, and HCO + . In this first scenario, then, 
most of the Arp 220 luminosity is produced in the 
ER, which has been modeled as a thin disk by 
SYB97, and as a warped disk by Eckart & Downes 
(2001). Significantly, if this model is correct, 
a spatially extended star burst, responsible 
for the bulk of the far-infrared luminosity 
is inferred. 

Our best fit to the continuum using 
model Si assumptions is presented in Figure 
5a, with derived physical parameters listed in Ta- 
ble 1. In all models, uniform densities throughout 
the different components are assumed for simplic- 
ity. In Table 1, X t is the wavelength for which 
the nucleus becomes optically thin (r = 1); ow- 
ing to the high opacities involved, the inferred 
emission is rather insensitive to the spectral index 
/?. Thus we have given values of the physical pa- 
rameters for ft = 1.5 and ft — 2. The dust mass 
has been derived by assuming a mass-opacity co- 
efficient of 12 cm 2 /gr (Hildebrand 1983). The 
parameters that have been allowed to vary in our 
models of the ER are the dust temperature Td, the 
diameter d (within the range 500-800 pc), and the 
dust opacity; 0 is fixed to 2 to ensure negligible 
emission at millimeter wavelengths (see S99). 

Model S 2 : In the second scenario (hereafter S 2 ), 
the emission from the nuclei at 25 /zm is assumed 
to be attenuated by foreground local dust with 
T"abs(24.5/zm) = 1.2, a value which was cho- 
sen to be compatible with the silicate ab- 
sorption observed in S99. The nuclei account 
for the required flux at 24.5 and 1300 /zm with 
Id % 106 K, significantly warmer than the 
~ 85 K temperature in Si and, as before, 
the emission in the submillimeter is opti- 
cally thick (see Table 1). In this scenario as well, 
however, the flux at 60-100 microns is again 
underestimated, and an emitting ER must also be 
involved to account for it. Nevertheless, the lu- 
minosity from the warm nuclei in S 2 is enough to 
account for the observed total luminosity from Arp 
220, so that the ER merely re-radiates the emis- 
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sion from the nuclei and no extended starburst is 
then needed to provide the bulk of the luminosity. 

Figure 5b shows our best fit for S 2 - The unique 
parameters that have been allowed to vary in our 
models of the EE are the diameter d and the dust 
opacity; 0 is again fixed to 2, and the dust temper- 
ature Td (shown in the insert panel of Figure 5b) 
has been computed from the requirement that the 
heating balance the cooling throughout the source. 
We have assumed spherical symmetry with the nu- 
cleus, the primary heating source, located at 
the center of the ER, which is divided into a set of 
spherical shells to account for the variation of Td 
with the distance to the nucleus. Despite the good 
fit to the continuum in Fig. 5b, S 2 implic- 
itly supposes a lack of spherical symmetry 
or some dumpiness, because the derived ra- 
dial opacity of the ER at 24.5 pm is 11.3, whereas 
the adopted opacity of the absorbing shell in front 
of the nuclei is r a bs(24.5/xm) = 1.2. Therefore 
these results should be considered approximate, 
but suggestive. The intrinsic geometry that un- 
derlies S 2 implies that the total luminosity, which 
coincides with the luminosity of the nuclei, is 
lower than the value inferred in Si , where spherical 
symmetry and uniformity is assumed for each 
component (Table 1). 

Table 1 shows that 0 — 2 yields a dust mass 
of ~ 10 8 M© for the nucleus, and that Si gives a 
mass 1.4 times higher than S 2 . Since SYB97 infer 
a dynamical mass of 6 — 8 x 10 9 M© enclosed in 
the inner 250 pc radius, and this region contains 
the nuclei and most of the inner disk (the ER), S 2 
with 0 = 1.5 is favoured in our models provided 
that the gas-to-dust mass ratio is not lower than 
the standard value of sslOO. 

It is worth noting that these models may be 
applied to the source as a whole, as implicitly as- 
sumed above, or alternatively to each one of an 
ensemble of N c smaller clouds of radius R c that 
do not spatially overlap along the line of sight. 
The value of N c x R 2 C determines the absolute 
scale, and the radial opacity and temper- 
ature distribution of each cloud as a func- 
tion of the normalized radial coordinate, i? c , 
determine the continuum shape. Identical 
results are found as long as the above pa- 
rameters remain constant. Furthermore, both 
alternatives give identical total masses, but dif- 
fer in the inferred mean density, which scales as 


y/N~ c . N c = 1 gives the lowest mean density 
< n(H 2 >, which is listed in Table 1. Typical 
values of a fewxlO 4 cm -3 are derived for the nu- 
clei, accounting for the emission from molecules 
with high dipole moment such as HCN, HCO + , 
HC 3 N and CN (see also SYB97). For the ER we 
obtain n(H 2 ) < 10 3 cm -3 ; since HCN and HCO + 
appear to show extended low-level emission (Rad- 
ford et al. 1991b), it is suggested that the actual 
density is higher than this lower limit or that the 
gas is clumpy. 

Although both scenarios S\ and S 2 reproduce 
the continuum emission from Arp 220 over the 
25-1300 pm interval and support the constraints 
on the nuclear sizes derived from the available 
high angular resolution continuum measurements, 
the dynamical masses inferred from CO millime- 
ter line observations favour S 2 with 0 = 1.5 for 
the nuclei . Moreover , as we discuss in sections 
5.3 and 6.1, the observed line ab sorption /emission 
also favours model S 2 . We thus adopt scenario S 2 
for the detailed modeling and analysis of the line 
emission and absorption. 

4.1. Extinction 

In both scenarios, the high brightness and com- 
pactness of the nuclei imply extreme continuum 
optical depth, corresponding to Ay ~ 10 4 mag. 
This conclusion is in strong contrast with the much 
more moderate values derived from infrared and 
radio hydrogen recombination lines (Genzel et al. 
1998; Anantharamaiah et al. 2000). The high ex- 
tinction derived here is the direct result of the 
measured 1.3 mm continuum flux from the nuclei, 
210 mJy (Carico et al. 1992; Sakamoto et al. 1999), 
and the observed upper limit of the correspond- 
ing source size, ~ 0'/4 (Sakamoto et al. 1999). 
These values imply r\ 3mm x Td{ K) ~ 40, which 
shows that even assuming unexpectedly high av- 
erage dust temperatures (e.g. Td = 200 K) and 
0=1 the dust emission is still optically thick 
even at 200 pm. On the other hand, the radio 
recombination lines observed by Anantharamaiah 
et al. (2000) are not affected by dust obscuration, 
although their predicted fluxes and the derived ex- 
tinction may be somewhat model dependent. These 
very different extinction values may be under- 
stood if we assume that the observed H recom- 
bination lines, tracing primarily star formation, 
are formed in the outermost regions of the nu- 
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Table 1 

Physical parameters derived from the continuum emission 


Scenario 

Component * 

'Tabs 

(24.5 /im) 

d 

(pc) 

P 

T d 

(K) 

(/im) 

M d 

(M©) 

< n(H 2 ) > a 
(cm- 3 ) 

L 

(L 0 ) 

Si 

Nucleus 

0 

145 

1.5 

85 

640 

5.1 x 10 7 

6.5 x 10 4 

4.8 x 10 11 





2.0 

85 

770 

1.3 x 10 8 

1.6 x 10 5 

4.8 x 10 11 

Si 

ER 

0 

630 

2.0 

50 

110 

5.2 x 10 7 

8.1 x 10 2 

8.3 x 10 11 

s 2 

Nucleus 

1.2 

145 

1.5 

106 

520 

3.6 x 10 7 

4.6 x 10 4 

9.8 x 10 11 





2.0 

106 

650 

9.3 x 10 7 

1.2 x 10 5 

9.8 x 10 u 

s 2 

ER 

0 

660 

2.0 

40-90 b 

82 

3.9 x 10 7 

5.3 x 10 2 

— 


a Assuming a single source instead of an ensemble of clouds 

b The temperature is calculated throughout the ER from the balance between heating and cooling with the 
nucleus as the primary heating source 



Kest (MHl) 

Fig. 3. — Line assignments in the vecinity of the 
153.2 /im feature. The grey line shows the adopted 
continuum level. 


clei, while a buried central energy source , respon- 
sible for the heating of dust in the innermost re- 
gions of the nuclei and a significant fraction of 
the galactic luminosity , is weak in recombination 
lines . If weak in recombination lines , the buried 
energy source is presumably weak in PAH features 
and PDR lines as well, consistent with the strong 
[C II] deficit in Arp 220 (see also section 5). Dust- 
bounded ionized regions, in which most of the Ly- 
man continuum from nuclear starbursts or AGN is 
absorbed by dust rather than by gas, may explain 
these properties of Arp 220 , as was proposed by 
Luhman et al. (2003). We further argue that the 
Lyman continuum luminosities derived from re- 
combination lines do not empirically rule out the 
possibility that an AGN accounts for more than 
~ 50% of the bolometric luminosity of Arp 220, 
because of the high dispersion of Lboi/TL y c val- 
ues shown by both starburst galaxies and AGN , 
the range of L L yc values derived from different 
tracers in Arp 220, AND THE UNCERTAINTIES 
IN THE ASSUMED EXTINCTION LAW 
AND THE DERIVED EXTINCTION (Gen- 
zel et al. 1998). Our derived N( H 2 )~ 10 25 cm“ 2 is 
high enough to obscure a source of high 5-10 keV 
luminosity from one or both nuclei in Arp 220, so 
that a hidden AGN is allowed despite the rela- 
tively weak X-ray luminosity observed in Arp 220 
(Clements et al. 2002). Haas et al. (2001) have 
also argued that a hidden AGN powers much of 
the luminosity of Arp 220 on the basis of the ob- 
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Fig. 5. — Fits of the continuum emission from 
Arp 220 in the 25-1300 fim range for scenarios (a) 
Si and (b) S 2 . Solid line shows ISO-LWS spec- 
trum of Arp 220, and filled triangles show the 60 
and 100 fim IRAS fluxes for comparison. Filled 
squares show the fluxes measured by Soifer et al. 
(1999, 24.5 /im); Eales, Wynn- Williams, & Dun- 
can (1989, 450 /mi); Rigopoulou, Lawrence, & 
Rowan-Robinson (1996, 350, 800, and 1100 /mi); 
and Sakamoto et al. (1999, 1300 /im). Dotted lines 
indicate the computed contributions from the nu- 
clei and the ER, whereas the solid grey line show 
the expected total flux. The insert panel shows the 
calculated dust temperatures in function of the ra- 
dial angular distance for S 2 (see text). 
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served submillimeter continuum excess relative to 
the 7.7 fim PAH flux. 

5. Models for molecules and atoms 
5.1. Comparison with Sgr B2 

The comparison of the spectrum of Arp 220 
with that of some well-studied galactic sources 
provides important clues about the regions where 
the observed lines are formed, while emphasizing 
the unique features that characterize the extra- 
galactic source. In this sense, Sgr B2 (component 
M) is an ideal comparison source, just because it 
shares common properties with Arp 220 de- 
spite the obvious differences in spatial scale 
(and indeed possibly in nature). Figure 6 shows 
the continuum- normalized spectra of Sgr B2 (M) 
(kindly provided by J.R. Goicoechea) and Arp 220. 
Sgr B2 harbors newly born OB stars, ultracom- 
pact H II regions, and hot cores, enshrined in a 
dusty envelope which is heated by shocks and by 
radiation, and which radiates a high IR luminos- 
ity (see GRC04 and references therein). The en- 
velope has the highest extinction in the direction 
of the N- and M-condensations, with optically 
thick emission in the far-infrared up to ~200 //m, 
and with foreground absorption lines of OH, HoO, 
CH, and [0 I] observed with the ISO/LWS grat- 
ing (GRC04, Fig. 6). Fabry-Perot observations 
of Sgr B2 have also allowed the detection of other 
molecular species like NH, NH 3 , NH 2 , HD, H 3 0 + , 
and C 3 , as well as high excitation lines of OH and 
H 2 O (GRC04 and references therein). As shown 
above (section 4), the extinction toward the nuclei 
of Arp 220 is also very high, although the dust is 
significantly warmer than in the Sgr B2 envelope. 
The comparison between both sources is mean- 
ingful (at least as a first approximation, and from 
the point of view of the radiative transfer) as long 
as the nuclear region of Arp 220 can be consid- 
ered an ensemble of continuum-thick molecular- 
rich clouds such as Sgr B2. For an ensem- 
ble of Sgr B2-like cloud, since both the continuum 
and the line absorption scale with the number of 
clouds, the continuum- normalized spectrum is the 
same as that of one individual cloud. This re- 
sult applies even if the lines in the ensemble are 
broadened relative to the one-cloud emission due 
to cloud-to-cloud velocity dispersion and rotation, 
provided that the lines remain unresolved with the 
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Fig. 6.— Continuum-normalized spectra of Sgr B2 (M) and Arp 220. The main carriers of some line features 
are indicated. 
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grating resolution. Thus the differences between 
the two spectra in Fig. 6 reveal real differences in 
excitation and/or column densities. 

The high-excitation lines of OH and H 2 O are 
much stronger in Arp 220 than in Sgr B2 (M) (see 
Fig. 6). In particular, the OH n 3 / 2 J = 9/2 — 7/2 
65 pm line, with strong absorption in Arp 220, is 
not detected in the grating spectrum of Sgr B2, 
and the OH II 3 / 2 J = 7/2 — 5/2 84 pm line is 
also much weaker in Sgr B2 (M). This strongly 
indicates the presence in Arp 220 of a high ex- 
citation region with relatively high OH column 
densities. Fabry-Perot spectral resolution allowed 
Goicoechea & Cernicharo (2002) to detect high ex- 
citation OH lines in Sgr B2, and showed that they 
are pumped through absorption of far-infrared 
photons. In Arp 220 also, these lines appear 
to be pumped by the strong infrared radia- 
tion flux in the neighbourhood of the nuclei 
(section 5.2). Toward Sgr B2 (M), Goicoechea 
& Cernicharo (2002) derived AT(OH)« 2 x 10 16 
cm -2 , and we may expect significantly higher col- 
umn densities toward the nuclei of Arp 220. The 
peculiarity of Arp 220 is also revealed by the rela- 
tively strong absorptions in the NH 3 and NH lines. 
These species have also been detected toward Sgr 
B2 (M) with Fabry-Perot spectral resolution (Cec- 
carelli et al. 2002, GRC04), but Fig. 6 indicates 
much higher column densities in Arp 220, at least 
toward the nuclei. 

GRC04 found that the OH 11* / 2 J = 3/2 - 1/2 
163 pm line, pumped through absorption of pho- 
tons in the OH 11^2 -n 3 / 2 J — 3/2 -3/2 53.3 pm 
line, shows emission over a large region associated 
with Sgr B2. In Arp 220 the line is strong, suggest- 
ing significant widespread emission, i.e. from the 
ER. It is also worth noting that, although the [C 
II] line is not detected in the grating spectrum of 
Sgr B2 (M), Fabry-Perot observations allowed its 
detection (Vastel et al. 2002, GRC04), with a flux 
of « 1.5 x 10“ 17 W cm -2 for the component 
observed in emission. This value is, within 
a factor of 2, similar to strengths in the 
surrounding region where the continuum is 
however much weaker (GRC04). Therefore, 
in addition to effects of self-absorption and 
absorption of the continuum by C + in fore- 
ground excitation clouds (Vastel et al. 2002, 
GRC04), the low [C II] /FIR ratio at Sgr 
B2 (M) is due to a strong increase in the 


FIR emission that is not accompanied by 
a corresponding rise of PDR line emission. 
Dust-bounded ionized regions (Luhman et 
al. 2003), together with extinction effects in 
the far- infrared, might account for the lack 
of PDR line emission associated with this 
additional infrared component. This is sim- 
ilar to our hypothesis for Arp 220’s nuclei 
(section 5.6). The [C II] emission in Arp 
220 is expected to arise from PDRs in the 
ER, where the bulk of the PAH emission 
is found (Soifer et al. 2002). As shown in 
Fig. 5, the ER dominates the observed far- 
infrared emission, but since the main heat- 
ing sources are, according to model S 2 , the 
nuclei, the ER is mainly heated via absorp- 
tion of infrared rather than UV photons , its 
intrinsic luminosity is relatively low, and 
hence the [C II] line remains weak (section 
5.6). 

Finally, we stress that Arp 220 and Sgr B2 (M) 
show similar absorptions in the ground-state lines 
of OH and H 2 0, as well as in the CH line at 149 
p m. Fabry-Perot observations of the low-lying OH 
and H 2 0 lines toward Sgr B2 (M) (Cernicharo et 
al. 1997, GRC04), indicate that most of these ab- 
sorptions are produced in diffuse low-excitation 
clouds located along the line of sight to Sgr B2 
but not physically associated with it. The simi- 
lar absorptions found in Arp 220 strongly suggest 
that a diffuse medium is also present there. In 
fact we have found that the combination model 
of the nuclei and the ER that reproduces reason- 
ably the high excitation OH and H 2 0 lines, fails to 
explain the strong absorptions in the lowest-lying 
OH and H 2 0 lines. The presence of an absorbing 
diffuse component is strengthened by the detection 
of CH at 149 /im, because most of the correspond- 
ing absorption toward Sgr B2, with strikingly sim- 
ilar strength, has been associated with translucent 
clouds (Stacey, Lugten, & Genzel 1987, GRC04). 
This component could also have an important ef- 
fect on the observed [0 I] 63 pm line: it is detected 
in absorption toward both sources and Baluteau 
et al. (1997) found that most of the absorption to- 
ward Sgr B2 could be produced by foreground gas 
in diffuse clouds. We will call this diffuse absorb- 
ing component “the halo”; the presence of such 
a halo in Arp 220 finds independent support in 
the outermost disk observed in the CO millimeter 
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lines by SYB97. 

5.2. Outline of the models 

The dust models described in section 4 set up 
the basis for the molecular calculations. These 
are carried out with the method described in 
Gonzalez-Alfonso & Cernicharo (1997, 1999), 
which computes the statistical equilibrium pop- 
ulations of a given molecule by assuming spherical 
symmetry and line broadening caused by micro- 
turbulence and/or radial velocity gradients. In the 
present calculations we have assumed, for simplic- 
ity, pure microturbulent line broadening, but some 
tests showed that the inclusion of a radial velocity 
gradient hardly modified the results. Rotational 
motion, which is present around the nuclear re- 
gion of Arp 220 (e.g. SYB97, Sakamoto et al. 
1999; Downes & Solomon 1998), is not included. 
Nevertheless, we note that for an arbitrary-shaped 
source rotating with uniform angular velocity 
around any axis, any two arbitrary points within 
the source have no Doppler shift along their join- 
ing segment, and therefore the statistical equilib- 
rium populations and emergent fluxes are identical 
to those computed by ignoring rotation. The in- 
clusion of uniform rotation would only broaden the 
expected emergent line profiles, for which our ISO 
spectrum provides no information. In the more 
general case of non-uniform rotation, our models 
may be considered a first-order approximation. 
The assumption of spherical symmetry and our 
neglect of steep velocity gradients that may result 
from cloud-to-cloud velocity dispersion MAY be 
considered more critical. 

Our non-local code accounts for radiative trap- 
ping in the molecular lines, collisional excitation, 
and excitation through absorption of photons 
emitted by dust. The dust parameters derived 
above for Si and S 2 are used in the calculations 
for molecules. As in the case of the continuum 
models, we have modeled the nucleus and the ER 
separately. This is required, in S 2 , by the rela- 
tively low dust opacities of the ER toward the 
nucleus as compared with the radial opacity of 
the ER (see above), and involves inevitably an 
additional uncertainty. 

Since the dust in the nucleus is very optically 
thick throughout the ISO wavelength range, model 
results are only sensitive to the molecular column 
densities in the most external parts of the nucleus. 


It is therefore assumed that only in the external- 
most regions of the nucleus, where the infrared 
lines are formed, the molecular abundances are 
different from zero (see Goicoechea & Cernicharo 
2002, for the case of Sgr B2). Since dust and 
molecules are assumed to be coexistent, extinction 
effects within the nucleus are implicitly taken into 
account; they place important constraints on the 
molecular abundances. We have adopted a molec- 
ular shell thickness of 2 x 10 18 cm, which for a 
mean n(H 2 )= 4.6 x 10 4 cm -3 (in S 2 ) corresponds 
to Ay ~ 50 mag and r(50/um) « 0.3. For lines 
around 50 ^m, the contribution to the absorption 
by molecules located deeper into the nucleus was 
checked in some tests to be relatively weak, due to 
both dust and molecular line optical depth effects. 
Foreground extinction in S 2 was also taken into 
account. In the models for the ERs, which have 
much lower continuum opacities, we have assumed 
that dust and molecules have uniform abundance 
ratio throughout the whole region. The presence 
of the central nucleus is included in the calcula- 
tion of the statistical equilibrium populations, but 
ignored in the calculation of the emergent fluxes 
to avoid accounting for it twice. 

Toward the nucleus, absorption of continuum 
radiation determines the excitation of OH and 
H 2 O; the radiative rates are much higher than 
the collisional ones even for molecular shock con- 
ditions (n(H 2 )= 5 x 10 6 cm' 3 and T*= 300 K). 
The rate coefficients of Offer, van Hemert, & 
van Dishoeck (1994) and Green, Maluendes, & 
McLean (1993) were used to check the collisional 
excitation of OH and H 2 O, respectively. If molec- 
ular shocks were present, only the absorption 
of the lowest-lying lines would be significantly af- 
fected. Since we use a halo to match these 
lines, we cannot distinguish between shock and 
non-shock conditions (i.e., the line ratios are not 
sensitive to n(H 2 ) and T* within plausible values). 
For simplicity, physical parameters for the halo are 
derived by assuming non-shock conditions for the 
nucleus. Concerning the ER, widespread shock 
conditions are not applicable because they would 
involve strong emission in the H 2 0 2 i 2 — loi and 
3 o 3 - 2i2 lines and in most OH lines (hard to can- 
cel by any halo). Our molecular data are therefore 
only sensitive to the radial molecular column den- 
sity N and the microturbulent velocity dispersion 
so that only these two computational parame- 
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ters are required to define a model of a given com- 
ponent. 

Similar to the dust models, models for 
molecules may be applied to the source as a whole, 
or alternatively to each one of an ensemble of 
smaller clouds that do not spatially overlap along 
the line of sight. Besides the scaling relationships 
pointed out above, the molecular column density 
must remain the same to obtain identical results 
when varying N c - On the other hand, once the 
number of continuum sources N c is fixed (e.g., 
N c — 1), the same line fluxes are obtained if 
both N and o v are divided by the same factor / c , 
and the resulting fluxes are then multiplied by / c . 
The latter reflects the approximate equivalence 
between one absorbing cloud with column density 
N and velocity dispersion <r v , and f c clouds, with 
parameters Nff c and cr v /f c , which overlap on the 
sky but not in the line-of-sight velocity space 1 . 

Variations in N and a v have different effects on 
the line absorptions. If a v decreases the absorp- 
tions are weaker for optically thick lines, but due 
to the increase of line opacities, the high-energy 
levels become more populated and the above weak- 
ness is more pronounced for the low-lying lines. 
On the other hand, the increase of N has little 
effect on very optically thick lines (most of them 
low-lying lines), but a larger effect on those with 
moderate opacities. 

We have generated a grid of models for the nu- 
cleus and the ER of both Si and S 2 , by varying the 
above free parameters N and o v . In each model, 
the molecular shell is divided into a set of sub- 
shells in order to account for the spatial varia- 
tion of the excitation temperatures of the lines. 
First we searched for the nucleus 4- ER combina- 
tion model that best matches the OH and H 2 O 
non-ground-state lines, with the same value of 
o v for both species. As pointed out above, the 
relatively deep absorptions of the OH and H 2 0 
ground-state lines could not be fitted satisfactorily 
by any model, and a halo component was added 
to match these lines. The halo was assumed to be 
a purely absorbing shell; although the equilibrium 
populations were computed in spherical symme- 
try assuming a size of three times that of the ER, 


Hn fact the models could be defined in terms of the two 
independent variables N/a v and N x / c ; nevertheless we 
will use the variables N and a v with / c = l. 


limb emission (i.e. emission for impact parame- 
ters that do not cross the continuum source) was 
ignored in the calculation of the fluxes emergent 
from the halo. 

Once the model for OH and H 2 O is determined, 
models for the other detected species are per- 
formed, also keeping fixed the value of a v derived 
above for each component. Figure 7 compares the 
derived model spectrum and the observed one, and 
Table 2 lists the inferred parameters. The next 
sections are devoted to explaining the details of 
these calculations and results. 

5.3. OH and H 2 0 

Our best model fits for the high excitation lines 
of OH and H 2 0 come from model S 2 , and the re- 
sults presented here will focus on this scenario. 
The main difference between Si and S 2 consists in 
the higher column densities and/or broader line 
widths required by Si to reproduce the lines, ow- 
ing to the fact that Td is significantly lower in 
Si . Models for the nucleus with broad line widths 
( 1 a v > 60 kms -1 , see below), however, predict 
strong absorptions in some OH and H 2 0 lines, 
such as the OH J = 5/2 - 3/2 98 /im and 
the H 2 0 3i3 - 2 0 2 1 38 /im lines, which are not 
observed. Therefore S 2 , which still requires high 
column densities, is favoured by our models. 

Our best model fits for the high excitation lines 
involve column densities of 2 — 6 x 10 17 cm -2 for 
both OH and H 2 0, and a v = 50 - 30 kms -1 , 
towards the nucleus. The model for the nucleus 
reproduces nearly the whole absorption in the OH 
IT 1/2 7/2 — 5/2 and ni / 2 — n 3 / 2 5/2 — 5/2 lines, 
most of the OH n 3 / 2 9/2 — 7/2, and significant ab- 
sorption in the other lines but by far too weak 
in the ground-state lines. It also reproduces 
the full absorptions in the H 2 0 4 22 - 3 i 3 , 4 32 - 3 2 i , 
3 3 o - 2 2 i , 3 3 i - 2 20j and 4 23 - 3 i2 lines, and signif- 
icant absorption in the 3 2 i - 2i 2 and others. The 
somewhat low value derived for a v is required to 
keep the absorptions weak in some low excitation 
lines which are marginally or not detected. Since 
low a v implies low velocity coverage for absorp- 
tion of the continuum, the column densities that 
are needed to explain the absolute values of the 
absorptions in high excitation lines are relatively 
high. We stress that this value of o v must be in- 
terpreted as a strong lower limit on the linewidths 
that would be observed with high enough spectral 
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Continuum-subtracted spectrum (10 



50 60 70 8 



Fig- 7. — Continuum-subtracted spectrum of Arp 220, compared with the result of the model (in grey). The 
line features that contribute more to the model are identified (see text). 
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Table 2 

Column densities toward Arp 220 a 


Species Nucleus ER Halo 


OH 

2.0 

X 

10 17 

2.0 

X 

10 17 

1.8 

X 

10 16 

h 2 o 

2.0 

X 

10 17 

3.0 

X 

10 16 

1.5 

X 

10 16 

CH 



- 



- 

2.0 

X 

10 15 

NH 

1.3 

X 

10 16 



- 

1.8 

X 

10 15 

nh 3 

2.8 

X 

10 16 



- 

4.0 

X 

10 15 

C+ 

2.0 

X 

10 19 

2.0 

X 

10 19 



- 

0 

5.0 

X 

10 19 

5.0 

X 

10 19 

2.7 

X 

10 18 


a Units are cm 2 


resolution, because rotation and cloud-to-cloud ve- 
locity dispersion would broaden the observed lines. 
The CO kinematic models of Arp 220 by SYB97 
found a generic value of a v = 90 km s” 1 ; as the au- 
thors discuss this value should be considered the 
joint effect of the local linewidth and the cloud- 
to-cloud velocity dispersion over a scale of ~100 
pc. Our o v is the local linewidth involved in the 
calculation of opacities and directly related to the 
molecular excitation, and thus the kinematic value 
of SYB97 must be considered here an upper limit. 

It is worth noting that the above column densi- 
ties are derived by forcing the OH and H 2 0 high 
excitation lines to arise in the same region. A 
slight improvement to the fit of the H 2 0 lines 
is obtained with even lower cr v , 25 kms” 1 , and 
N( H 2 0)~ 10 18 cm -2 . This may indicate that 
HoO and OH toward the nucleus do not arise in 
the same regions. We will adopt in the following 
the nucleus model with a v — 50 kms* 1 , corre- 
sponding to the spectrum shown in Fig. 7 and the 
parameters given in Table 2, and estimate an un- 
certainty of a factor of ~3 on the derived column 
densities. 

Once the model for the nucleus is matched, a 
search of the combination ER+halo that better ac- 
counts for the remaining flux in non ground-state 
lines is searched for. The ER mainly accounts 
for the OH Hi/ 2 3/2 — 1/2 163 [i m line, which 
is uniquely observed in emission (but predicted in 
absorption in the nuclei), for more than half of the 
absorption observed in the OH n 3 / 2 7/2 -5/2 line, 
for reemission in the fl 2 / 2 5/2 - 3/2, and for sig- 


nificant absorption in the three ground-state OH 
lines. We estimate for the ER iV(OH)~ 2 x 10 17 
cm” 2 , with o v — 50 kms” 1 throughout most of 
the ER and o v — 90 kms” 1 just around the nu- 
cleus. This higher value of a v was required to ob- 
tain significant reemission in the non-detected OH 
n 3/2 7/2 — 5/2 line; since here geometrical effects 
may be important, this result should be consid- 
ered with caution. For H 2 0 we obtain a signifi- 
cantly lower iV(H 2 0)~ 3 x 10 16 cm” 2 , giving sig- 
nificant absorption in the 3 22 — 3 i 2 , 2 20 — In, and 
221 — 1 10 lines, and some reemission in the 3o3 — 2 12 
line. Since the halo also yields some absorption in 
these H 2 0 lines but much deeper absorption in the 
ground-state 2 22 - l 0 i one (see below), the rela- 
tive H 2 0 column density in these components is 
not well determined. 

In the halo, the values of N( OH) and a v were 
determined by fitting the missing absorption in the 
three OH ground-state lines. A value of a v rela- 
tively low, 15-20 kms” 1 , and iV(OH)~ 2 x 10 16 
cm" 2 , were found to reasonably fit the fluxes of 
the 79 and 119 gm lines, though the flux of the 
53 ji m line is somewhat underestimated (Fig. 7). 
Higher values of a v would predict too much ab- 
sorption in the already saturated 119 /xm line. For 
H 2 0 we find A(H 2 0)~ 1.5 x 10 16 cm" 2 . Owing 
to the strong radiation field from the nucleus and 
ER, H 2 0 in the halo is still significantly excited, 
thus yielding also some absorption in the 2 2 o — In, 
221 — lio and 3 0 3 - 2 i2 lines. 

Despite the generally satisfactory fit obtained 
for the OH and H 2 0 lines, the shoulder of the 
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OH 119 fim line, presumably produced by the 
18 0H II3/2 5/2 - 3/2 line, is not reproduced. The 
models assume 16 OH/ 18 OH=500; however, the 
isotopic abundance ratio in Arp 220 required to 
reasonably fit the 120 /im shoulder is as low as 
16 OH/ 18 OH^50. This value cannot be ruled out 
if 14 N is efficiently converted into 18 0 in nuclear 
processing of high mass stars, and then efficiently 
ejected to the interstellar medium through stel- 
lar winds and/or supernovae (Henkel & Mauers- 
berger 1993). Also, there is compelling evidence 
for isotopic ratios of 150-200 in starburst regions 
of nearby galaxies (Henkel & Mauersberger 1993). 
Due to the low spectral resolution of the spectrum, 
and the possibility that the feature is contami- 
nated by other species, we do not attempt to place 
useful constraints on this ratio. Nevertheless, we 
conclude that a very low 16 OH/ 18 OH abundance 
ratio likely applies to Arp 220, perhaps indicating 
an advanced stage starburst (Henkel & Mauers- 
berger 1993). 

5.4. CH 

We assume that the CH line we observe at 149.2 
fim arises in the halo, based on the results ob- 
tained toward the Galactic Center (see GRC04). 
We derive iV(CH)« 2 x 10 15 cm' 2 by fitting the 
feature. Unlike the case of Sgr B2, however, the 
far-infrared emission from the nuclear region of 
Arp 220 is so strong that significant excitation of 
the ( N , J)=(l,3/2) rotational level occurs, and 
our models for the halo predict some contribu- 
tion by CH to the absorption at 181 fim (Fig. 7). 
The latter is uncertain, however, because the fea- 
ture at 181 fim. could be contaminated by H3O+ 
(GRC04) and/or by stronger absorption of H^O, 
whose abundance relative to H^O could be en- 
hanced relative to the assumed value of 1/500. 
On the other hand, there is a wing-like feature 
at 118.5 fim, observed in both the forward and re- 
verse scans, which could be caused by the doublet 
CH (3,7/2)— (2.5/2). If so, and since the excitation 
of this line requires a relatively strong radiation 
field, there would be CH in the nuclei that would 
account for about 1/3 of the absorption at 149.2 
fim, and the CH column density in the halo would 
be 2/3 of the quoted value. 


5.5. NH and NH 3 

Assuming that the 153.2 fim absorption fea- 
ture is caused entirely by NH, we have obtained 
7V(NH)~ 10 16 cm -2 toward the nucleus. The 
model reproduces the marginal absorption feature 
at 76.8 fim, attributable to the Nj = 4j - 3 
lines. For the above column density, the 2 2 — li 
and 2 3 — I2 lines at 153.2 fim become saturated 
and the associated feature is not completely repro- 
duced, so that we have added an additional halo 
component with N(NH)« 2 x 10 15 cm” 2 (Table 
2). Nevertheless, this model still underestimates 
the absorption at 102 fim, strongly suggesting the 
contribution from other species (section 3). 

Within a given K— ladder, the excitation of the 
NH 3 non-metastable levels is determined by 
absorption of far-infrared continuum pho- 
tons, while the metastable levels in K = 
2,3,... are pumped through collisions (see 
e.g. Ceccarelli et al. 2002, for the case of 
Sgr B2). In the model of Fig. 7 we have as- 
sumed an average density of n(H 2 )= 4.6 x 10 4 
cm -3 and 7* = 100 K, but we have checked 
that the model results are insensitive to 
the adopted T* because of the blending of 
lines from different K— ladders to each spec- 
tral feature in our spectrum (Fig. 1). We 
obtain iV(NH 3 )~ 3 x 10 16 cm” 2 toward the nu- 
cleus to fit the absorptions at 125 and 127 fim 
(Fig. 7). Besides the absorption in the lines 
showed in Fig. 4, the model predicts significant 
absorption at «100 fim and «101.6 fim, caused 
by {J,K) = (5 ,K) - (4 ,K) lines. The model, 
however, fails to explain the strong absorption at 
166 fim, and therefore an halo component with 
jV(NH 3 )~ 4 x 10 15 cm -2 has been added to the 
global model of Fig. 7. Nevertheless, the halo 
components of NH and NH 3 should be con- 
sidered uncertain , because variations in the 
background continuum associated with each 
component could in principle account for 
the missing flux in the lines. 

Finally, we have explored the possibility that 
the NH 2 radical contributes to the spectrum at 
some wavelengths. The expected strongest ab- 
sorption from NH 2 is found at ssl59.5 fim, caused 
by the strongest components of the 3 i 3 - 2 02 or- 
tho line (the hyperfine structure was neglected in 
these calculations, but the split of the levels due 


17 

Appendix C 



to the unpaired electronic spin of 1/2 was taken 
into account). At this wavelength, a marginal ab- 
sorption feature may be attributed to NH 2 , and 
is approximately fit with a model for the nucleus 
where N( NH 2 ) 10 15 cm 2 . We have included 

this model of NH 2 in Fig. 7 to show that the ex- 
pected absorption in other lines, like the 3 22 - 2 U 
one at 105 /zm, do not conflict with the observa- 
tions, and we conclude that iV(NH 2 )< 2 x 10 15 
cm -2 . 

5.6. C II and O I 

A crucial test of our model is whether it can 
reproduce the [C II] 157.7 /im emission and [0 I] 
63.3 /im absorption lines. Among ULIRGs, Arp 
220 shows one of the most extreme [C II] deficits 
(Fcii/Ffir - 2x 10 -4 , Luhman et al. 2003). The 
[C II] line is formed within Ay <2 mag from the 
surfaces of PDRs (e.g. Wolfire, Tielens, & Hol- 
lenbach 1990), where the UV field from nearby 
high mass stars, or from the average galactic field 
has not been significantly attenuated. In this re- 
gion, photodissociation maintains most of the gas 
in atomic or singly ionized form, but some radicals, 
like OH and NH, find their maximum abundances 
there (Sternberg & Dalgarno 1995). In particu- 
lar, OH is expected to be an excellent molecular 
tracer of PDRs’ surfaces, given that its abundance 
is rather low in UV-shielded quiescent molecular 
clouds. Its abundance relative to H nuclei within 
the C + region of dense PDRs is expected to ap- 
proach the value of ~ 3 x 10“ 6 (Sternberg & Dal- 
garno 1995). In fact, Goicoechea & Cernicharo 
(2002) have found an abundance of » 2 x 10 -6 in 
Sgr B2, and it could be as high as 5 x 10~ 6 around 
the galactic center (Genzel et al. 1985). We have 
estimated X(OH)~ 1 - 3 x 10 -6 toward the nu- 
cleus of Arp 220. 

On the above grounds, and adopting a gas 
phase carbon abundance of 1.4 x 10“ 4 (Savage & 
Sembach 1996), we assumed N(C~ h )/N(OH)=100 
in Arp 220 (Table 2), and we computed the ex- 
pected [C II] line emission by assuming excitation 
through collisions with H (Tielens & Hollenbach 
1985). The collisional rates were taken from Lau- 
nay & Roueff (1977a), and H densities of 5 x 10 4 
cm“ 3 and 10 3 cm -3 were assumed for the nucleus 
and ER, respectively. Since the critical density is 
3 x 10 3 cm -3 (Kaufman et al. 1999), results for 
the nucleus are not critically dependent on the as- 


sumed density. They are also insensitive to the 
assumed temperature as long as it is higher than 
~100 K (Wolfire et al. 1990). Thus our results are 
only sensitive to the assumed density in the 
ER, and to the assumed C + abundance. 

The result of this calculation has also been 
added to the overall model of Fig. 7 ( T* = 500 
K has been assumed), providing a reasonable fit 
to the [C II] line emission. The contribution from 
the nucleus is only 14% of that from the ER be- 
cause of the low volume of the nuclear emitting 
region. The expected line flux from the nucleus 
could be lower if absorption of the underly- 
ing continuum by low-excitation C + in the 
halo, ignored in this model, occurs as ob- 
served in Sgr B2 (M) (GRC04). The bulk of the 
line emission arises from the ER. The situation 
toward Arp 220 resembles what is found in Sgr 
B2 (GRC04), where the line is emitted mainly 
from an extended region around condensa- 
tions N and M, while the strong FIR source 
itself is not associated with corresponding 
observable [C II] line emission. 

The [O I] 63.3 /zi n line has been modeled by as- 
suming N ( O I)/iV(OH)=250 (Table 2); although 
the oxygen abundance is expected to be twice that 
of C+ in the atomic region, it is expected that fur- 
ther atomic oxygen exists deeper into the clouds 
(Sternberg & Dalgarno 1995). For this reason, a 
wide range of excitation temperatures is expected 
for the [O I] line. We have just fitted a single 
“effective” kinetic temperature and assumed also 
collisions with atomic H. The same densities as as- 
sumed above for C + excitation are used for O I, 
and the collisional rates are taken from Launay & 
Roueff (1977b). 

Our calculations show that absorption in the 
[O I] 63.3 /zm line is obtained, both toward the 
nucleus and the ER, with an effective T* = 200 
K, but the line is still too weak to account for the 
observed feature. Given the high O I abundances 
that are expected in diffuse clouds (Baluteau et 
al. 1997), we have added to the model a halo com- 
ponent with N(0 I)= 3 x 10 18 cm” 2 (Table 2). 
This model accounts for the observed absorption 
at 63.3 /zm (Fig. 7). 

On the basis of the low extinction derived from 
infrared and radio H recombination lines and the 
high optical depth derived from our dust models 
(section 4.1), and the assumption that [C II] emis- 
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sion is expected to suffer extinction similar to that 
of the recombination line emission , we infer that 
extinction does not play a role in the / C II] deficit 
and favor the explanation proposed by Luhman 
et al. (2003, see also section 5.1) that non-PDR 
far-infrared emission is responsible for the [C II] 
deficit in Arp 220. Nevertheless, our continuum 
models by themselves (section 4) cannot rule out 
the possible role of far-infrared extinction on the 
measured line fluxes: the derived high nuclear far- 
infrared opacities indicate that only the most ex- 
ternal regions of the nuclei, where the OH and 
H2O lines are formed, and the ER, are able to 
contribute to the [C II] line emission. Our model 
fits cannot discern whether the intrinsic nuclear 
[C II] emission is negligible or rather is obscured 
by dust. In either case, we have shown that the 
[C II] line is well reproduced by assuming that C + 
and OH are coexistent, 

6. Discussion 

6.1. The Extended Region (“ER”) 

Our models support the widespread presence 
of PDRs in the ER. Both the high OH and C 
II column densities indicate that the UV field 
from newly formed stars have a profound effect on 
the chemistry in the ER. Significant contributions 
from shocks can be neglected, as was pointed out 
in section 5.2. The H 2 0-to-0H abundance ratio is 
significantly lower than 1, probably indicating en- 
hanced H 2 0 photodissociation. Both the models 
for dust and the [C II] 158 /zm line emission indi- 
cate rather moderate densities (nn ~ 10 3 cm" 3 ). 

Our models also indicate that the bulk of the 
[C II] line emission arises in the ER. It is there- 
fore likely that star formation is responsible for 
this emission. This result is strongly supported 
by the observations of Soifer et al. (2002), who 
found that the PAH emission is also spatially ex- 
tended. Assuming a “ normal ” [C II] /FIR ratio, 
i.e. F[c ii]/Ffir = 5x 10~ 3 (Stacey et al 1991 ) , 
the expected intrinsic FIR emission from the ER 
is ~ 3 x 10 10 Lq, i.e. ~3% of the total galactic 
luminosity. This estimate supports the scenario 
*5*2 that has been used to model the line emission. 
According to our models , the total FIR luminosity 
(due to absorption and re-emission of nuclear in- 
frared radiation ) from the ER is much higher than 
the intrinsic (PDR) emission. Our estimate of the 


intrinsic ER luminosity is somewhat lower than 
that derived by Soifer et al. (2002) and Spoon et 
al. (2004), w ho found a lower limit of ~ 7 x 10 10 
L 0 , but confirms qualitatively their results. The 
starburst luminosity in the ER seems to be similar 
to that of moderately bright infrared galaxies, like 
NGC 253 (Radovich, Kahanpaa, & Lemke 2001). 

The warm mass in the ER derived from our 
OH models is ~ 10 9 Mq, assuming X(OH)=10~ 6 
relative to H nuclei and the warm H 2 arises from 
the same volume as the OH. This estimate of mass 
may be considered an upper limit. From the H 2 ro- 
tational lines detected by ISO, Sturm et al. (1996) 
estimated ~ 3 x 10 9 M 0 of warm gas. Given the 
uncertain correction by extinction, it is likely that 
an important fraction of the H 2 emission arises in 
the ER. 

6.2. The halo 

The models also indicate the presence of a 
halo. Despite the uncertainties in the column 
densities of this component, they are typical of 
those found toward our galactic nucleus: iV(H 2 0, 
OH)~ 2 x 10 16 cm -2 have been also derived in 
the diffuse medium toward Sgr B2 (Cernicharo et 
al. 1997; Neufeld et al. 2003; Goicoechea & Cer- 
nicharo 2002), and the derived 7V(CH)« 2 x 10 15 
cm -2 is also similar to that found toward Sgr B2 
by GRC04 and Stacey et al. (1987). The column 
densities in the halo derived for NH and NH3 are 
uncertain because they are based on single lines, 
but they could also exist in a population of molec- 
ular clouds located far away from the nuclear re- 
gion. 

6.3. The nuclei 

The similarity of the OH column densities in 
the ER and the nucleus may suggest that, at 
least to some extent, we are observing the same 
widespread OH component, and that OH is more 
excited toward the nucleus because of the underly- 
ing stronger infrared continuum in that direction. 
On this ground a PDR origin of the observed OH 
would be favoured. However, the inferred high 
H 2 0 column density would be difficult to explain 
in this context. Although H 2 0 is expected to form 
efficiently in UV-shielded regions of dense PDRs 
(Sternberg & Dalgarno 1995), with total column 
densities similar to those of OH, there seems to be 
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no clear correlation between iV(OH) and 7V(H 2 0). 
In fact, the H 2 0 column density in the ER is sig- 
nificantly lower than that of OH. Also, detection 
of NH3 and, above all, of NH, seem to indicate an 
additional nuclear component. 

OH could also arise in C-shocks (Watson et al. 
1985), but it is unlikely that they dominate the 
OH absorptions because iV(H 2 0) would then be at 
least one order of magnitude higher than A (OH). 
However, an interesting possibility is that those 
C-shocks, or alternatively hot core regions, are 
combined with PDRs, i.e., H 2 0 produced there 
is subject to a strong UV field and thus to pho- 
todissociation. This process may be enhanced if 
there is hot gas emitting X-rays, like in supernova 
remnants where the OH abundance is expected to 
be > 10“ 6 (Wardle 1999), or from a nuclear AGN. 
In the nuclei of Arp 220, several compact 18 cm 
continuum sources indicate the presence of high lu- 
minous supernovae (Smith et al. 1998); however, 
the extended and external OH required to ex- 
plain the infrared data suggests a more widespread 
component. The diffuse OH megamaser emission 
found in Arp 220 should be related to it (Lonsdale 
et al. 1998). Soft-extended and hard-compact X- 
ray emission, detected around and from the nuclei 
(Clements et al. 2002), could be responsible for 
photodissociation of H 2 0 produced in shocks and 
hot cores, thus enhancing the OH abundance. In 
particular, hot cores are expected to exist widely 
in the nuclei, given the high dust temperatures 
and densities found there. The presence of J- 
shocks, where iV(OH) is expected to be higher 
than N(R 2 0) except for high enough preshock 
densities (Neufeld & Dalgarno 1989), cannot be 
disregarded. 

The high column densities obtained for NH and 
NH 3 seem to indicate that standard gas-phase 
PDR chemistry alone is not able to explain the full 
molecular data in Arp 220. Sternberg & Dalgarno 
(1995) predicted for a PDR an NH column density 
more than two orders of magnitud lower than that 
of OH, whereas we estimate 7V(OH)/A r (NH)~20 
in Arp 220. The enhancement of NH relative to 
OH may be more than one order of magnitude. 
In Sgr B2, GRC04 also found a somewhat high 
iV(NH) relative to OH, i.e. 7V(OH)/iV(NH)~30- 
100. It is interesting that the high NH abundance 
in diffuse clouds is a factor of ~ 40 higher than 
predicted by gas-phase chemical models (Meyer & 


Roth 1991). The latter has been used by Wagen- 
blast et al. (1993) and Crawford & Williams (1997) 
to argue for grain-surface production of NH. In 
principle, this process could also help to enhance 
the NH abundance in Arp 220, because the dust 
in the nuclei has been found to be warm so that 
grain mantles could efficiently evaporate, and also 
because in an enviroment with enhanced cosmic 
rays, as expected from a starburst and the conse- 
quent high rates of supernovae, or an AGN, the 
release of grain mantles to the gas phase by means 
of sputtering should be also enhanced. However, 
hydrogenation of NH should in principle continue 
until saturation, because indeed we observe NH3 
with a column density twice of NH. The issue 
that now arises is that, if hydrogenation gener- 
ally completes, then a very low NH-to-NH 3 abun- 
dance ratio would be expected, and if it does not 
complete like in the models of Wagenblast et al. 
(1993), the scenario fails because of the low rel- 
ative abundance found for NH 2 . In Sgr B2, for 
example, NH:NH 2 :NH 3 =1:10:100 (GRC04), and 
NH 2 is found to be fairly abundant (van Dishoeck 
et al. 1993). 

One possible solution for the low NH 2 abun- 
dance is that the already invoked enhancement of 
cosmic rays deeply affects the ion-molecule gas- 
phase chemistry. Federman, Weber, & Lambert 
(1996) have shown that the high NH abundance 
found in some diffuse galactic enviroments could 
be explained through cosmic ray ionization of 
atomic nitrogen, followed by hydrogen abstrac- 
tion reactions that form NH+ and dissociative re- 
combination that yields NH. The slight endother- 
micity of N + -hH 2 -4NH + -hH (Millar, Farquhar, & 
Willacy 1997) is not a problem here, given the high 
dust temperatures in the nuclei. If in Arp 220 the 
cosmic-ray ionization of N is enhanced, the above 
scheme could give rise to high NH abundances. 
Furthermore, both NH 2 and NH 3 would be much 
less abundant if photodissociation is important in 
those regions. NH 3 would be formed primarily 
in grain mantles through nearly complete hydro- 
genation, thus again keeping the NH 2 abundance 
low, and released to the gas phase in widespread 
hot core regions relatively shielded from UV fields. 
H 2 0 could also follow this last process. 

A chemistry deeply influenced by ion-neutral 
reactions have been also invoked by Aalto et al. 
(2002) to explain the high emission from HNC rel- 
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ative to HCN in Arp 220 and other luminous in- 
frared galaxies. Furthermore, Aalto et al. (2002) 
also found strong subthermal CN emission, indica- 
tive of gas at moderate densities and irradiated by 
UV fields (Rodriguez-Franco, Martin-Pintado, & 
Fuente 1998). In the most simple scenario, the 
molecular content of Arp 220 may thus be inter- 
preted in terms of hot cores submitted to strong 
UV and X-ray fields, where enhanced evaporation 
of grain mantles and ion-molecule chemistry in- 
duced by cosmic ray ionization are also deeply af- 
fecting the relative molecular abundances. 

7. Summary 

We have analyzed the ISO/LWS spectrum of 
Arp 220 using radiative transfer models applied to 
both the continuum and line emission. Our main 
results are: 

1. The continuum emission from 25 to 1300 /zm is 
well reproduced with a two-component model: (a) 
the nuclei, with effective size of 0"4 and dust tem- 
perature of 106 K, which accounts for essentially 
the whole flux at 25 fim and at millimeter and 
submillimeter wavelengths, and (b) an extended 
region (ER), whose effective size is 2" and which 
dominates the continuum emission from 60 to 250 
fi m. 

2. The extinction toward the nuclei is very high 
(Ay ~ 10 4 mag); the dust in the ER is heated 
through absorption of radiation emanating from 
the nuclei. 

3. The spectrum of Arp 220 shows molecular lines 
of OH, H 2 0, CH, NH, and NH 3 , as well as the 
atomic [O I] 63 /zm line in absorption and the 
[C II] 158 fim line in emission. The outermost 
regions of the nuclei, along with the ER, 
are traced by the lines observed in the far- 
infrared. 

4. The high excitation lines of OH and H 2 0 are 
pumped through absorption of photons emitted by 
dust. Column densities of iV(OH)= 2 - 6 x 10 17 
cm -2 and 7V(H 2 0)= 2-10x 10 17 cm" 2 are derived 
toward the nuclei. In the ER, N( OH)~ 2 x 10 17 
cm" 2 and iV(H 2 0)^ 3 x 10 16 cm" 2 . We found it 
necessary to invoke a third component, or halo, 
to match the low-lying lines of OH and H 2 0; 
this halo has column densities that are similar to 
those found toward the Galactic Center (iV^OH, 
H 2 0)~ 1.5 x 10 16 cm" 2 ). 

5. The CH line detected in the far-infrared spec- 


trum of Arp 220 is assumed to arise from the 
halo, and the inferred column density is V(CH)~ 
2 x 10 15 cm" 2 . This value is also similar to that 
found toward the Galactic Center. 

6. Models for NH and NH 3 indicate high column 
densities toward the nuclei, iV(NH)~ 1.5 x 10 16 
cm" 2 and iV(NHs)~ 3 x 10 16 cm" 2 . The upper 
limit found for the column density of NH 2 is much 
lower, iV(NH 2 )< 2 x 10 15 cm" 2 . 

7. The [C II] 158 fim line strength is approxi- 
mately reproduced by assuming that C + is 100 
times more abundant than OH. Our models pre- 
dict that the line arises mainly from the ER, and 
that non-PDR far-infrared emission , with possible 
extinction effects, is mostly responsible for the ob- 
served [C II] deficit in Arp 220 . The [O I] 63 /zm 
line is also matched with an abundance of 250 rel- 
ative to OH and absorption toward the nuclei, the 
ER, and the halo. 

8. PDR molecular chemistry plays a key role in 
the ER and contributes to the elevated OH abun- 
dance at the expense of Zf 2 O. Toward the nuclei, 
however, important contributions from hot cores, 
and possibly from shocks, is most plausible. The 
nitrogen chemistry, and in particular the high NH 
abundance, seems to be strongly influenced by ion- 
neutral reactions triggered by cosmic ray ioniza- 
tion. 
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ABSTRACT 

The prototypical Seyfert type 2 galaxy NGC 1068 has been observed with the 
Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory 
(ISO, ESA). This provides the first complete far-infrared spectrum from 43 to 
197/im. In addition to the 7 expected ionic fine structure emission lines, we also 
detect OH rotational lines at 79, 119 and 163/mi. Unlike any other galaxy yet 
observed, NGC 1068 shows all these lines in emission. To increase the information 
on the strong ionic emission lines powered either by the Seyfert nucleus and 
the starburst, we have included in our analysis measurements in the 2-35 ^m 
spectral range from the ISO Short Wavelength Spectrometer (SWS), as well as 
optical and ultraviolet lines from the literature. To model the observed line 
intensities we have considered three emission components: the AGN component, 
the starburst component originated from the circumnu clear ring at ~ 15 " from 
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the center and the nuclear star burst component. Following the guidelines of a 
recent work, the nuclear emission line spectrum from the AGN was modeled using 
two components photoionization models with constant densities of 2000 cm -3 and 
10 4 cm -3 , respectively, and emission regions extending from a distance from the 
center of ~ 20 pc to ~ 200 pc. We have adopted three different ionizing continua 
illuminating the AGN narrow line region: the first one has a large absorption dip 
around 4 Ryd; the second has the standard broken power-law shape generally 
used for NGC1068 and the third has a big blue bump peaking at ~ 10 Ryd. We 
find that the UV to mid-IR emission line spectrum originated from the AGN 
cannot be used to severely constrain the ionizing continuum that illuminates the 
NLR of NGC1068, contrary to previous results. While the presence of a big blue 
bump is ruled out, our modeled spectra do not require the presence of a large 
dip at 4 Ryd., but are consistent also with the canonical power-law type ionizing 
spectrum. 

We have used the Br 7 emission from the circumnuclear ring and that one 
originated from the nuclear region to normalize the starburst emission of the 
second and third component, respectively. The ionizing continua for the starburst 
models have been derived from starburst synthesis models. 

Three of the FIR lines of OH have been detected in NGC1068, and we present 
the first results of our modeling, which suggests the gas lies in small clouds, about 
O.lpc in size, with densities of < 10 4 cm~ 3 at the outer edges. 

Subject headings: galaxies: individual (NGC 1068) - galaxies: active - galax- 
ies: nuclei - galaxies: Seyfert - galaxies: emission lines - galaxies: starburst - 
infrared: galaxies. 


1. INTRODUCTION 

NGC 1068 is among the most studied active galaxies. It is known as the archetype 
Seyfert type 2 galaxy, it is nearby, with a measured redshift z=0.0038 (Huchra et al. 1999) 
(corresponding to a distance of D=15.2 Mpc, adopting H 0 =75 km s -1 Mpc -1 ), providing a 
scale of only ~ 74 pc/", bright (L IR = 2 x 10 n L© Bland-Hawthorn et al. 1997) and it has 
been extensively observed and studied in detail from X-rays to radio wavelengths. Around 


HSO is an ESA project with instruments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom) and with the participation of ISAS and NASA. 
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the active nucleus, NGC 1068 has a central nuclear star cluster extending ~0.6" (Thatte et 
al. 1997). Moreover, a stellar bar of 2.3kpc of size has been observed in the near-IR (Scoville 
et al. 1988; Thronson et al. 1989), which is surrounded by a circumnuclear starburst ring. 
Mapping revealed that the infrared emission in NGC 1068 was due to both the Seyfert nucleus 
(dominating the 10/wn emission) and to the star forming regions in the bright circumnuclear 
ring of ~ 3kpc in diameter (emitting most of the luminosity at A > 30 pm) (Telesco et al. 
1984). A Br'y imaging study (Davies, Sugai Ward 1998) reveals a similar morphology. 
It indicates that a short burst of star formation occurred throughout the circumnuclear 
ring (15-16 " in radius) perhaps within the last 4-40 Myr. CO interferometer observations 
revealed molecular gas very close to the nucleus (~0.2") suggesting the presence of ~ 10 8 M© 
within the central 25pc (Schinnerer et al. 2000). Recent high resolution H 2 line emission 
mapping indicates the presence of two main emission knots with a velocity difference of 140 
km/s, that, interpreted as quasi-keplerian velocity, would imply a central enclosed mass of 
1O 8 M 0 (Alloin et al. 2001). 

In this article, we present the first complete far- infrared spectrum from 43 to 197^/m 
showing both atomic and molecular emission lines (§2.). We model the far-IR continuum 
emission using gray body functions (§3.), we also model the composite UV- to far-IR atomic 
emission line spectrum, as derived by combining literature observations with our data, using 
photoionization models of both the active nucleus (§4.1) and the starburst components (§4.2- 
4.3). Moreover, montecarlo radiation transfer models have been used to account for the 
detected far-IR molecular lines (§5.). Our conclusions are then given in §6. 


2. OBSERVATIONS 

NGC 1068 was observed with the Long Wavelength Spectrometer (LWS) (Clegg et al. 
1996) on board the Infrared Space Observatory (ISO) (Kessler et al. 1996), as part of the 
Guaranteed Time Programme of the LWS instrument team. The full low resolution spectrum 
(43-197/rm) of NGC 1068 was collected during orbit 605 (July 13, 1997). Two on-source 
full scans (15,730 seconds of total integration time) and two off-source (6’ N) scans of the 
[CII]l58yum line (3,390 seconds of total integration time) were obtained. On- and off-source 
scans had the same integration time per spectral step. Because of the design of the LWS 
spectrometer, simultaneously with the 158//m data, a short spectral scan of equal sensitivity 
to the on-source spectrum was obtained at sparsely spaced wavelengths across the LWS 
range. 

The LWS beam is roughly independent of wavelength and equal to about 80 arcsec. 
The spectra were calibrated using Uranus, resulting in an absolute accuracy better than 
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30% (Swinyard et al. 1996). The data analysis has been done with ISAP 2 , starting from 
the auto-analysis results processed through the LWS Version 7-8 pipeline (July 1998). To 
be confident that newer versions of the pipeline and calibration files did not yield different 
results, we have compared our data with the results obtained using pipeline 10.1 (November 
2001) and we did not find significant differences on both the line fluxes and the continuum. 
To increase the signal-to-noise ratio, all the full grating scans taken on the on-source position 
have been coadded, as well as the two sets of data on the off-source position. No signal was 
detected in the off-source position. The emission line fluxes were measured with ISAP, 
which fits Gaussians on top of a local continuum. In every case the observed line width was 
consistent with the instrumental resolution of the grating, which is typically 1000 km/sec. 
The integrated line fluxes measured independently from data taken in the two scan directions 
agreed very well, to within 10%. LWS has a known tendency for fringing to produce spurious 
ripples in the continuum baseline (Swinyard et al. 1998). Fortunately, these fringing ripples 
are hardly noticeable in our LWS spectrum, presumably because the far-IE continuum is 
centrally concentrated towards the center of the LWS 80 arcsec beam. 

Besides the LWS observations, we also use the SWS observations presented by Lutz 
et al. (2000), to extend the wavelength and ionization-level coverage. Table 1 presents all 
the ISO line flux measurements, including those from the SWS. 


3. THERMAL CONTINUUM SPECTRUM 

The 50-200/im far-infrared spectral energy distribution (SED) of NGC1068 can be fit- 
ted by the combination of different gray body functions. Fig. 1 shows the observed LWS 
spectrum of NGC1068, together with a fit with the main gray body component at T=32K, 
a warmer component at T=65K and a cooler one at T=20K. The dust emissivity has been 
assumed to be proportional to the inverse square of the wavelength (i.e., the flux density is 
S\ — B\(T) x A _/9 with Ba(T) the Planck function at the wavelength A and temperature T 
and /?=2). We have adopted this steep dust emissivity law because it fits better than flatter 
laws the 12-200 //m SED of Seyfert and starburst galaxies, as shown in Spinoglio, Andreani 
& Malkan (2002). 


2 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are Centre d’Etude Spatiale des Rayonnements (France), 
Institute d’Astrophysique Spatiale (France), Infrared Processing and Analysis Center (United States), Max- 
Planck-Insitut fur Extraterrestrische Physisk (Germany), Rutherford Appleton Laboratories United King- 
dom) and the Space Research Organization, Netherlands. 
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Because of the large (3) number of components used, we do not claim that the proposed 
fit is unique, however it fits the 50-200/xm spectrum fairly well (reduced x 2 =0.40, corre- 
sponding to a probability of P>98%) and it is derived from a physically consistent view of 
the various dust thermal components in Seyfert galaxies (see e.g. Spinoglio, Andreani & 
Malkan 2002). 


4. THE FINE STRUCTURE LINES 

To be able to better constrain the modeling of the line emission of NGC1068, we have 
combined our far-infrared fine structure line measurements (Table 1) with ultraviolet, optical 
and infrared spectroscopic data from the literature (Kriss et al. 1992; Marconi et al. 1996; 
Thompson 1996; Lutz et al. 2000). The complete emission line spectrum of NGC1068 from 
the ultraviolet to the far-IR includes several low-ionization lines that are primarily produced 
outside the narrow line region (NLR) of the active nucleus, as well as intermediate ionization 
lines that originate from both starburst and AGN emission. For this reason, we find that 
no single model satisfactorily explains all the observed emission lines. We identify three 
components: 

- an AGN component (the NLR), exciting the high ionization lines and contributing 
little to the low-to-intermediate ionization lines; 

- a starburst component in the circumnuclear ring of the galaxy (e.g. Davies, Sugai &: 
Ward 1998) that produces part of the emission of the intermediate ionization lines; 

- another starburst component that originates from the nucleus of NGC1068. 

We will also see that a photo-dissociation region (PDR) (e.g. Kaufman et al. 1999), in 
the interstellar medium of the galaxy, has also to be added to fully model the low ionization 
and neutral forbidden lines. 

The contributions of the starburst and PDR emission can be estimated from the CII/- 
OI/OIII line ratio diagram, as used by Malkan et al. (2002). However, for this galaxy we 
can obtain a more direct measurement of the starburst/HII region components, using the 
high spatial resolution available from ground-based near-IR Fabry-Perot imaging and grating 
spectroscopy. As shown in §4.2 and §4.3, we have taken the Br 7 fluxes originated in the ring- 
like star forming complex and from the nuclear region of NGC1068, respectively from Davies, 
Sugai & Ward (1998) and from Oliva &: Moorwood (1990). 
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In the following of this section, we will examine separately the three components that 
are producing the composite fine structure emission line spectrum of NGC1068, for which 
we propose three different computations in §4.4. We will leave for §4.5 the estimate of the 
model parameters of the PDR emission produced from the ISM of the galaxy. 


4.1. Modeling the AGN 

The first predictions based on photoionization models of the mid to far-infrared spec- 
trum originated from the Narrow Line Regions (NLR) of active galaxies have been done by 
Spinoglio & Malkan (1992), well before the observations could be collected. The specific 
photoionization modeling of only the high ionization emission lines unambiguously produced 
by the Seyfert nucleus of NGC1068 has been done previously by Alexander et al. (2000). 
Following their results, we have simulated their models, although using a different photoion- 
ization code, CLOUDY (Version 94.00 Ferland 2000), and then we have varied the shape 
of the ionizing continuum. As already mentioned in the previous section, our goal was not 
only to test if the Alexander et al. (2000)’s results were unique, but also to add together all 
different photoionized components originated from different emission regions to model the 
complete emission line spectrum of NGC1068. We have assumed two component models for 
the AGN, with the same parameters as in Alexander et al. (2000): the first component has 
a constant Hydrogen density of 10 4 cm~ 3 , an ionization parameter U=0.1, extends from ~ 
20 to ~ 80 pc from the center, it has a covering factor c=0.45 and a filling factor of 10 -4 
with a radial dependence of the form r~ 2 ; the second component has a density of 2 x 10 3 
cm -3 , an ionization parameter U=0.01, extends from ~ 140 to ~ 190 pc from the center, 
it has a covering factor c=0.29 and a filling factor of 10 -3 without any radial dependence. 
We have assumed also the ’’low Oxygen” abundances adopted by Alexander et al. (2000), 
to be able to compare in the easiest way our results with theirs. Tables 2 and 3 report the 
results of the three AGN models we have chosen to fit the NGC1068 spectrum, together 
with the observed line fluxes: the line fluxes are given for each of the two components A 
and B, which are treated as independent, and the total flux for each model is simply the 
sum of the fluxes of the two components. The intrinsic nuclear spectrum of NGC 1068 has 
been inferred by Pier et al. (1994), we have modified this in the 1-100 Ryd region to model 
both the through derived by Alexander et al. (2000) and a big blue bump, as expected by the 
thermal emission of an accretion disk around the central black hole. Specifically, model A 
assumes the same ionizing spectrum that was derived from the best fit model of Alexander 
et al. (2000), i.e. with a deep through at 4 Rydberg (log/ = -27.4, -29.0, -27.4, -28.2 at 
2, 4, 8 and 16 Ryd, respectively). Model B assumes the original ionizing spectrum derived 
from Pier et al. (1994). Model C assumes a Big Blue Bump spectrum on the top of the 
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Pier et al. (1994) ionizing continuum (log / = -25.8, -25.8, -25.8, -27.4 at 2, 4, 8 and 16 Ryd, 
respectively). The three AGN ionizing continua are plotted in Fig.2. The inner and outer 
radii of the emission regions of the two components, 20, 80, 140 and 190pc, correspond to 
angular distances of about 0.27, 1.1, 1.9 and 2.6 ", respectively. 

We can see from Tables 2 and 3 that only the AGN A and AGN B models , and not 
the AGN C model, represent fairly well most of the observed line fluxes, except for the low 
and intermediate ionization lines, that — as we expected — are also, and sometimes mainly, 
produced by starburst and PDR excitation (see Section 4.4). This first result rules out the 
presence of a big blue bump in the ionizing continuum of NGC1068. To be able to compare 
the modeled ultraviolet and optical lines with the observations, we also listed in Table 2 and 
3 their dereddened fluxes, assuming two values for the extinction: E b-v = 0.4 mag (Malkan 
&; Oke 1983) and E b-v — 0.2 mag (Marconi et al. 1996). 


4.2. Modeling the starburst ring 

NGC1068 is known to emit strong starburst emission from both its nucleus and the 
ring-like structure at 15-16 "from the nucleus, traced e.g. by the Br 7 emission (Davies, 
Sugai &; Ward 1998). To model this second component, we have considered how much of the 
Br 7 line was emitted inside the different apertures of the ISO spectroscopic observations. 
In the smallest ISO-SWS aperture used (14 " x 20 ") no Br 7 emission is detected from the 
ring. In the next aperture (14 " x 27 ") Davies, Sugai &: Ward (1998)’s regions I, J, K and L 
are included, for a total Br 7 intensity of ~ 10 -14 erg s -1 cm -2 . The largest SWS aperture 
(20 " x 33 "), in addition to the above regions, also contains regions H, E, D and C, for a 
total Br 7 intensity of ~ 2.2 x 10 -14 erg s _1 cm -2 . Finally, in the large LWS aperture (~ 
80 " 2 ), all the ring-like Br 7 emission is contained, ~ 5.9 x 10 -14 erg s _1 cm -2 . 

Assuming that applies ’’case B” recombination, we can scale the Br7 intensity into Ha 
intensity (from Table 4.4 of Osterbrock 1989): 


jBr y /jH a - JBt-, / jH a • 3 Br y / 3 H e ‘ 3hJ 3H 0 ’ jHpl 3 H a = 0.0097 

For the three increasing apertures (14 " x 27 ", 20 " x 33 " and 80" 2 respectively) we 
have derived an Ha intensity from the starburst ring of 10~ 12 , 2.3 x 10 -12 and 6.1 x 10 -12 
erg s _1 cm -2 , respectively. We have used these values to derive the modeled line intensities 
from the CLOUDY code, that was used adopting Ha as normalizing line. 

We have chosen three starburst synthesis models Starburst99 from Leitherer et al. 
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(1999), as input ionizing spectra for the CLOUDY photoionization code. These ioniniz- 
ing continua are shown in Fig. 3. For the models with instantaneous star formation law, 
we adopted a total mass of M = 10 6 M 0 , and a duration of 5 Myr, while for the models 
with continuous star formation law, we have chosen a star formation rate of 1 M© yr -1 . For 
both types of models we adopted a Salpeter IMF (a=2.35), a lower cut-off mass of 1 M©, 
an upper cut-off mass of 100 M©, abundances of Z=0.008 and nebular emission included. 

We report in Table 4 the line fluxes predicted using CLOUDY for the two instantaneous 
models (SBR A and B, where SBR stands for Starburst Ring) and the four continuous models 
(SBR C, D, E and F). The adopted parameters for the photoionization models are given in 
the table. 


4.3. Modeling the nuclear starburst 

We have computed also the expected emission that originates from the nuclear starburst 
present in NGC1068. To model this third component, we have used the Br 7 intensity of 13 
x 10~ 14 erg s _1 cm -2 , which was observed in a 6"x 6"aperture (Oliva & Moorwood 1990). 
As a first order approximation, we made the assumption that half of the Br 7 intensity 
observed in the nuclear region is originated from the nuclear starburst, while the other half 
is produced in the NLR. Then we used the two models F and E, described in the previous 
section. The predicted line fluxes are listed in Table 3, together with the model parameters. 


4.4. Adding the three components 

Summing up the line intensities of each one of the three components, the composite 
spectrum of NGC1068 can be derived and compared with the observed one. We have chosen 
three combinations to compute the composite models, each one with a different AGN model: 
1) the first one (that we name CM1, as Composite Model 1) with the AGN ionizing continuum 
as suggested by Alexander et al. (2000) (model AGN A), a low-ionization nuclear starburst 
(model SBN F) and an high-ionization starburst ring (model SBR E); 2) the second (CM2) 
with the original Pier et al. (1994) (model AGN B), the same nuclear starburst and starburst 
ring models as for CM1; 3) the third (CM3) with the hypothetical bump (model AGN C), 
and again the same starburst models. The results of these three composite models are given 
in Table 5, compared to the observed and dereddened values, assuming the two choices 
for the extinction (see §4.1). We also show the results of the three composite models in a 
graphical way in Fig. 4 and Fig. 5, where the modeled to the observed flux ratio is given 
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for each line and for the two values assumed of the extinction, respectively. 

Taking the higher value for the optical extinction (Fig.5), most of the lines (33/43) of the 
first model (CM1) and almost all the lines (37/43) of the second composite model (CM2) are 
predicted within a factor 3 from the observed/dereddened values. The third composite model 
(CM3) is much less successful, having 14 lines badly under- or over-predicted. Considered 
the fact that our composite models are made of three different physical components, aimed 
to reproduce more than forty emission lines, the agreement between the first two composite 
models and the data is satisfactory and, as already noted, excludes the possibility that a big 
blue bump is present in the ionizing spectrum of NGC1068. On the other hand we cannot 
confirm the presence of the dip suggested by Alexander et al. (2000). 


4.5. Modeling the PDR 

An accurate model for the PDR emission originating in the ISM of the NGC1068 galaxy 
is beyond the scope of the present article, however it is possible to check what are the main 
parameters of the PDR, i.e. far-ultraviolet field (Go) and Hydrogen density (n). To do this, 
we just subtract off the modeled intensities of the three components from the observed fluxes 
of [OI]63pm, [OI]145pm and [CII] 158pm, then we take the ratios of the these lines, namely 
[OI]63pm/[CII] 158pm and [OI]63pm/ 145pm and compare them with the Kaufman et al. 
(1999) PDR models. Following this simple procedure, we find that for the first composite 
model log Go = 2 and log n = 3.3, while for the other two composite models log Go = 3 and 
log n = 2.5. 


5. THE OH LINES 

In NGC 1068, we detect three of the OH rotational lines. Two of the lines are to 
the ground state 2 n 3 / 2 3/2 level: the in-ladder 119pm fundamental line from the 2 n 3 / 2 5/2 
level, and the cross ladder 79pm line from the lowest, 2 Hi/ 2 1/2 level. The third line is the 
lowest transition of the 2 ni/ 2 ladder: the 163pm line between the J=l/2 and J=3/2 levels. 
These three lines are all seen in emission in NGC 1068, in striking contrast with the OH 
emission we have measured from other bright infrared galaxies such as Arp 220 (Fischer et 
al. 1999), Mrk 231 (Harvey et al. 1999), NGC 253 (Bradford et al. 1999), M 82 (Colbert 
et al. 1999), in which the 119 micron fundamental is in absorption (the 79pm lines is seen 
both in emission and absorption; the 163pm line is always seen in emission). In addition 
to these three detections, our ISO LWS observations were able to obtain useful limits on 
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four of the other eleven FIR OH lines in the LWS band that arise between the lowest eight 
rotational levels. The lines are indicated in the overall LWS spectrum of NGC1068 in Fig. 1, 
are shown in detail in Fig. 6 and tabulated in Table 1, while a comparison between observed 
and modeled fluxes is given in Table 6. 

We used a montecarlo radiative transfer code to model the lines in NGC 1068, as part 
of a detailed study of all the OH lines observed by ISO in galaxies (Smith 2002; Smith et 
al. 2003). The code was developed by the Submillimeter Wave Astronomy Satellite (SWAS) 
mission (Ashby et al., private com.), and is a modification of the original Bernes code, but 
includes a treatment of continuum photons from dust mixed in with the gas - a particularly 
essential feature for OH, which is pumped in many cases by absorption of 35 fxm continuum. 
In addition the code corrects for some previous errors encountered at large optical depths, 
also an issue of importance for OH which has a very strong matrix element (1.67 Debye), 
and a fundamental line transition rate of about 0.1 s -1 - about 1000 times greater than CO 
rates. The code also includes an ability to handle a wide range of molecules besides rigid 
rotors. 

The montecarlo code is one-dimensional; so far this constraint does not seem to preclude 
a reasonable physical explanation. The code takes in input a series of concentric shells, 
each of which is specified as to size, gas and dust temperature, H 2 density, velocity and 
turbulent velocity width, and molecular abundance relative to H 2 . The model allows the 
wavelength dependence of the dust emissivity to be independently set. The SWAS version 
of the montecarlo code calculates the populations of the molecular levels in each shell, and 
the output of this calculation, iterated until it converges, is then fed into a radiative transfer 
code that calculates the line profiles as seen by an external observer looking at the cloud. 

The montecarlo output by itself is not constrained to be self-consistent; for example, 
the input parameters need not conserve luminosity between shells. To obtain this self- 
consistency, which we determined was essential for a satisfactory model of the continuum 
flux, we used the DUSTY code (Ivezic k Elitzur 1997) to model the continuum and generate 
a set of shell parameters that provided this consistency. We then iterated the DUSTY 
model with the montecarlo line outputs. Although we do not claim this technique gives a 
unique solution for the cloud structure, it does give a canonical model consistent with the 
observations. 

We take as a starting assumption that the clouds of gas and dust producing the OH 
emission are heated from the inside, rather than dominated by external heating; the most 
successful models implied the cooler volumes were substantially larger than the hotter vol- 
umes. To produce consistent emission, especially in the 119 (xm fundamental line, we find 
that the clouds are small and dense. The best fit agreement has clouds of dimensions of 
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~0.2 pc in radius, and densities of < 10 4 cm -3 at the outer edges, increasing towards the 
center with a power law behavior of ( R/Rq ) 1-25 . The temperature in the outer shells is about 
22K, increasing inwards approximately with a power law dependence (R/Ro) 0 68 . Our best 
fit model has a total column density of H 2 in each cloud is 1.5 x 10 24 cm -2 , and a relative 
abundance of OH in the model is 10 -7 ; note that we take OH to be absent in the hot, inner 
portions of the cloud (T > 300K) . The strong lines of OH are therefore very optically thick, 
and precise radiative transfer calculations like the present modeling are essential. 

As can be seen from Table 6, the montecarlo code predicts the line flux ratios to the 
119/xm feature to within a factor of 40% of the observed values. The limits to the observed 
fluxes in all the weaker lines are also consistent with this modeling. At a distance of 15.2 Mpc, 
the total number of clouds needed to produce the observed absolute flux in the lines and 
the continuum is ~ 3 x 10 7 . As for the continuum, the montecarlo code also successfully 
reproduces the continuum emission seen from NGC 1068, to within a factor of 2-3 in absolute 
flux density across the entire LWS spectrum, under the simple assumption that the dust is 
1% of the gas by mass, and that it has the same temperature as the gas everywhere in the 
cloud. Although these constraints are not required by the code, they provide a good fit and 
are adequate for our general purposes here. The total mass in such an ensemble of clouds 
is approximately 9 x 10 9 M©. This amount of mass is considerably larger than the value of 
~ 10 8 M© estimated from CO interferometry in the core, but includes all of the galaxy in the 
ISO LWS beam of about 80", a much larger region than the circumnuclear CO observations. 
Papadopoulos & Seaquist (1999a) find that the mass in 1 'might be ~ 6 x 10 9 M©, or perhaps 
even twice that - so, in good agreement with our values. Our modeling of the large beam 
ISO flux averages is unable to discriminate between nuclear clouds and non-nuclear clouds. 
Moreover, Papadopoulos & Seaquist (1999b) also find a gas/dust ratio up to 300 in the 
central 40 ", much more than the standard value of 100 that was adopted here. 

The montecarlo code automatically calculates the velocity profile of each of the OH 
lines. We found that for a wide range of parameters the lines had both emission and ab- 
sorption components, the result of cold gas in the line-of-sight to the warm center absorb- 
ing, while the cooler, outer shells emitted. Our observations of the 119/rni, feature in the 
galaxy NGC 253 (Bradford et al. 1999), for which we have ISO-LWS Fabry- Perot velocity 
resolved measurements, reveal this dual character. The detection we report here, of signifi- 
cant emission in NGC 1068, does not preclude the possibility that the actual line emission 
is substantially larger, but has been offset by some absorption at the line center by other 
components. Furthermore the model also suggests that the absorptions we typically detect 
at low LWS spectral resolution in other ISO galaxies may in fact show, at higher resolutions, 
a combination of emission plus absorption features. 
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By way of a self-consistent check, we have also used the montecarlo code to model, inde- 
pendently, the emission and absorption from this series of clouds in the two neutral oxygen 
[01] lines at 63/xm and 145/xm. The results are satisfactorily consistent with the observa- 
tions. The [01] 63/im and 145/im lines from these warm compact clouds are, respectively, 
estimated at 1.2 x 10 -12 erg s -1 cm ~ 2 and 2.6 x 10 -14 erg s -1 cm~ 2 , which correspond to 
only a few percent of the total observed emission (see Table 1). Obviously since the [01] 
lines have their major contributions coming from other kinds of regions - photoionization, 
PDRs and shocks in particular - the significant conclusions to be drawn here are that the 
model for the OH lines is consistent with all the other observations. 

Our conclusion is that NGC 1068 has a large component of relatively small, dense clouds 
of molecular gas and dust, heated internally. We cannot say from these data where these 
clouds are located in the large ISO beam, or what if any relationship they may have with 
the nuclear regions except that they do not appear to require either pumping or radiative 
effects from the nuclear region or its torus. A further analysis of the OH lines, in the context 
of other galaxies measured by ISO, is in preparation (Smith et al, 2003). 


6. CONCLUSIONS 

The main results of this article can be summarized as follows: 

• The complete far-infrared (50-200^m) spectrum of NGC1068 has been observed for the 
first time. 

• The 50-200/zm continuum can be interpreted as thermal dust emission mainly due to 
three components at temperatures of 20, 32 and 65K, assuming a steep (/?=2) dust 
emissivity law. 

• The far-infrared ISO-LWS spectrum has been complemented with the mid-infrared 
data of ISO-SWS and with shorter wavelength (UV, optical and near-IR) literature 
data to assemble a composite atomic spectrum as complete as possible with the aim 
of modeling the different line emission components at work. This approach has been 
necessary especially because of the poor spatial resolution of the past space infrared 
spectrometers, which were not able to spatially separate the emission components. The 
lines have been interpreted as arising from three physically distinct components: the 
AGN component and two starburst components, the first one nuclear and the second 
one located in the ring at a radius of 15-16 " from the nucleus. Three composite 
models have been computed that vary mainly for their AGN component: the first one 
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has the ionizing continuum as derived from Alexander et al. (2000), showing a deep 
through at energies of a few Rydberg; the second has the monotonically decreasing 
ionizing continuum given by Pier et al. (1994) and the third has a big blue bump. Two 
values of the visual extinction ( Eb-v = 0.2 and 0.4) have been adopted to correct the 
optical and ultraviolet line fluxes. If we assume the higher value for the reddening, the 
agreement between the two composite models that do not show the presence of a big 
blue bump and the observed emission line spectrum is rather satisfactory, taking into 
account either the simplicity of the photoionization models chosen to avoid dealing 
with too many free parameters and the large number of lines which are originated in 
different physical regimes. 

• The radiative transfer modeling of the three lines of OH detected in NGC1068 suggests 
that the emission can be originated in small gas clouds, about O.lpc in size, with 
densities of < 10 4 cm -3 at the outer edges. 
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A(^m) 


Fig. 1. The ISO-LWS spectrum of NGC1068. The long-dashed line, which is fitting the 
data (reduced x 2 =0.40, corresponding to a probability of P>98%) represents the sum of the 
three different grey-body functions: dotted line: main component at T=32K; dotted-dashed 
line: cold component at T=20K; short-dashed line: warm component at T=65K. 
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E(Ryd) 


Fig. 2. — The AGN ionizing continua used as input for the photoionization models of 
NGC1068. The three continua differ in the frequency region between l>E/j y d>100, while 
outside this region the Pier et al. (1994) spectrum was adopted. The solid line shows the 
continuum derived from Alexander et al. (2000); the dashed line shows a simple power law 
interpolation; the dotted line shows the presence of the predicted ” big blue bump” . 


Appendix C 



- 18 - 



Fig. 3 — The starburst ring ionizing continua used as input for the photoionization models 
of NGC1068. The three continua are taken from Leitherer et al. (1999). 
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0.1 1 10 100 

transition X(/xm) 


Fig. 4.— The comparison of the composite models with the observations is shown as the 
ratio of modeled to observed flux ratio for each line, assuming a reddening of E(B-V)=0.2. 
Panels from top to bottom: model CM1, CM2 and CM3. The short dashed lines represent 
flux ratios within a factor 3 either ways. All line whose ratios exceed this factor are labeled. 
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Fig. 5. — Same as the previous figure, but with E(B-V)=0.4. 
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Fig. 6. — The observed emission lines of the OH molecule: from top to bottom: the 79^tm, 
the 119^m and the 163^tm lines. 
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Table 1. Measured line fluxes from the LWS and SWS grating spectra, with lcr 

uncertainties. 


Line 

A 

Flux 

Aperture 

reference 


(/im) 

(10 — 13 erg s~ l cm -2 ) 

("2) 


[Si IX] 3 P 2 -> 3 Pi 

2.584 

3.0 

14 x 20 

1 

[Mg VIII] 2 P 3/2 -» 2 P 1/2 

3.028 

11. ± 1.1 

14 x 20 

1 

[Si IX] 3 P, -f 3 P 0 

3.936 

5.0 ± 0.6 

14 x 20 

1 

[Mg IV] 2 P 1/2 -4 2 P 3/2 

4.487 

7.6 ± 1.5 

14 x 20 

1 

[Ar VI] 2 P 3/2 -4 2 P ]/2 

4.529 

15. ± 3. 

14 x 20 

1 

[Fe II] a4 P 9/2 -4 a6 Dg /2 

5.340 

5.0 

14 x 20 

1 

[Mg VII] 3 P 2 -4 3 P, 

5.503 

13. 

14 x 20 

1 

[Mg V] 3 Pi -4 3 P 2 

5.610 

18. ± 2. 

14 x 20 

1 

[Ar II] 2 P 1/2 -4 2 P 3/2 

6.985 

13. 

14 x 20 

1 

[Na III] 2 P 1/2 -»• 2 P 3/2 

7.318 

5.8 

14 x 20 

1 

[Ne VI] 2 P 3/2 -4 2 P 1/2 

7.652 

110. ± 11. 

14 x 20 

1 

[Fe VII] 3 P 4 -4 3 P 3 

7.815 

3.0 

14 x 20 

1 

[Ar V] 3 P 2 -4 3 Pi 

7.902 

< 12. 

14 x 20 

1 

[Na VI] 3 P 2 -4 3 Pj 

8.611 

< 16. 

14 x 20 

1 

[Ar III] 3 Pi -4 3 P 2 

8.991 

23.0 ± 3.3 

14 x 20 

1 

[Fe VII] 3 F 3 -*• 3 F 2 

9.527 

4.0 

14 x 20 

1 

[S IV] 2 P 3/2 -4 2 P 1/2 

10.510 

58. ± 6. 

14 x 20 

1 

[Ne II] 2 P 3/2 -4 2 P 1/2 

12.813 

70. 

14 x 27 

1 

[Ar V] 3 Pi -4 3 P 0 

13.102 

< 16. 

14 x 27 

1 

[Ne V] 3 P 2 -*• 3 Pj 

14.322 

97. ± 9.7 

14 x 27 

1 

[Ne III] 3 Pi -4 3 P 2 

15.555 

160. ± 32. 

14 x 27 

1 

[Fe II] a *Fi/i — ► a *Fg / 2 

17.936 

< 10. 

14 x 27 

1 

[S III] 3 P 2 -4 3 Pi 

18.713 

40. 

14 x 27 

1 

[Ne V] 3 Pi -4 3 P 0 

24.317 

70. ± 7. 

14 x 27 

1 

[O IV] 2 P 3/2 -4 2 Pj /2 

25.890 

190. ± 20. 

14 x 27 

1 

[Fe II] “ 6 D 7/2 -4 a 6 Dg / 2 

25.988 

8. 

14 x 27 

1 

[S III] 3 Pi -4 3 P 0 

33.481 

55. 

20 x 33 

1 

[Si II] 2 P 3/2 -4 2 Pi /2 

34.814 

91. 

20 x 33 

1 

[Ne III] 3 P 0 -4 3 Pi 

36.013 

18. 

20 x 33 

1 

[O III] 3 Pi -4 3 Pl 

51.814 

114. ± 3. 

80 

2 

[N III] 2 P 3/2 -4 2 P V2 

57.317 

51.4 ± 2.5 

80 

2 

[O I] 3 Pl -4 3 Pi 

63.184 

156. ± 1. 

80 

2 

[O III] 3 Pi -4 3 P 0 

88.356 

111. ± 1. 

80 

2 

[N II] 3 Pi -4 3 Pi 

121.897 

30.5 ± 1.1 

80 

2 

[O I] 3 P 0 -4 3 Pi 

145.525 

11.9 ± 0.4 

80 

2 

[C II] 2 P 3/2 -4 2 Pi /2 

157.741 

216. ± 1. 

80 

2 

OH 2 ni /2 5/2- 2 n 3/2 3/2 

34.60/34.63 

< 3. 

20 x 33 

2 

oh 2 ni /2 i/2- 2 ni /2 3/2 

79.11/79.18 

14.4 ± 1.5 

80 

2 

OH 2 n 3/2 5/2- 2 n 3/2 3/2 

119.23/119.44 

11.9 db 1.2 

80 

2 

OH 2 IIi /2 3/2- 2 ni /2 l/2 

163.12/163.40 

7.42 ± 0.65 

80 

2 


Note. — (1): from Lutz et al. (2000) and, where errors are available, Alexander et al. (2000); 
(2): this work 
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Table 2. Comparison of observed line fluxes with AGN model predictions 


Line id. A 

(Mm) 

Observed/D 1 /!} 2 

Flux (10- 13 erg s _1 cm" 
AGN A model 3 

Component A Component B 

- 2 ) 

AGN B model 4 

Component A Component B 

O VI A .1032+.1037 

37.4/4334./402. 

91.8+69.5 

3.0+2. 6 

14.5+11.1 

0.4+0.4 

(Lya)n A .1215 

101.8/3562. /602 

1014. 

1412. 

733. 

1327. 

N IV] A .1487 

5.1/103./22.9 

25. 

9.0 

23.6 

9.7 

(CIV)„ A .1549 

39.7/790./ 177. 

618. 

156. 

351. 

161. 

Hell A .1640 

21.4/426./95.5 

139. 

122. 

169. 

119. 

[Ne V] A .3426 

15.7/95./38.7 

103. 

34. 

58.5 

11.6 

[Ne III] A .3869+. 3968 

19. 2/97-/43.2 

10.4+3.1 

53.+16. 

17.3+5.2 

79.8+24. 

Hell A .4686 

6.1/27.6/13. 

19.2 

17.5 

23.9 

17.2 

[O III] A .4959+.5007 

256./964./496 

113.+39. 

424.+ 147 

91.5+31.6 

486. +168. 

[Si VI] A 1.96 

8.0/9.2/8.6 

9.3 

4.0 

9.0 

5.0 

[Si VII] A 2.48 

8.3 

5.1 

0.4 

6.8 

0.8 

[Si IX] A 2.584 

3.0 

0.7 


1.0 


[Mg VIII] A 3.028 

11. 

3.0 


4.6 


[Si IX] A 3.936 

5.4 

1.2 


1.8 


[Mg IV] A 4.487 

7.6 

4.5 

11.1 

6.7 

10. 

[Ar VI] A 4.529 

15. 

16.2 

14.5 

18.7 

12.7 

[Mg VII] A 5.503 

13. 

8.2 

0.5 

8.8 

0.3 

[Mg V] A 5.610 

18. 

13.4 

12. 

9.3 

6.1 

[Ar II] A 6.985 

13. 

0.1 

1.3 



[Na III] A 7.318 

5.8 


0.55 

0.1 

0.9 

[Ne VI] A 7.652 

110. 

111. 

17.5 

65.8 

5.3 

[Fe VII] A 7.815 

3.0 

3.8 

3.7 

3.2 

2.1 

[Ar V] A 7.902 

< 12. 

1.8 

3.5 

3.2 

4.0 

[Na VI] A 8.611 

< 16. 

1.2 

0.3 

0.8 

0.1 

[Ar III] + [Mg VII] A 8.991 

25. 

0.9+10. 

11.1+0.6 

1.2+11.0 

8.9+0.4 

[Fe VII] A 9.527 

4.0 

15.6 

15.8 

13.3 

8.9 

[S IV] A 10.510 

58. 

25.6 

86. 

29.1 

145. 

[Ne II] A 12.813 

70. 


7.3 


0.9 

[Ar V] A 13.102 

< 16. 

2.0 

5.1 

3.3 

6.1 

[Ne V] A 14.322 

97. 

132. 

98. 

102. 

45.9 

[Ne III] A 15.555 

160. 

3.4 

47. 

7.1 

67.9 

[S III] A 18.713 

40. 

4.5 

40. 

2.7 

28.3 

[Ne V] A 24.317 

70. 

58. 

77. 

42.3 

35.4 

[O IV] A 25.890 

190. 

38.5 

104. 

58.2 

105. 

[S III] A 33.481 

55. 

1.0 

22. 

0.6 

15.7 

[Si II] A 34.814 

91. 


1.5 


0.8 

[Ne III] A 36.013 

18. 

0.2 

4.0 

0.5 

5.8 

[O III] A 51.814 

110. 

1.8 

42. 

1.7 

48. 

[N III] A 57.317 

51. 

0.7 

14.5 

0.6 

13. 

[O I] A 63.184 

156. 




... 

[O III] A 88.356 

110. 

0.2 

10.2 

0.2 

11.6 

[N II] A 121.897 

30. 


0.2 



[O I] A 145.525 

12. 





[C II] A 157.741 

220. 


0.3 

... 



^ereddened line flux, assuming E b-v — 0.4 
2 Dereddened line flux, assuming Eq-v =0.2 

3 AGN A parameters: component A: Log U=-l., Log n=4, internal radius ~ 19 pc, external radius ~ 80 pc, ionizing 
spectrum from Alexander et al. (2000): component B: Log U=>2., Log n=3.3, internal radius ~ 140 pc, external radius 
~ 190 pc, ionizing spectrum from Alexander et al. (2000). 

4 AGN B parameters: same as AGN A models, but with the ionizing spectrum from Pier et al. (1994) 
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Table 3. Comparison of observed line fluxes with AGN and SBN model predictions 


Line id.A 
(*un ) 

Observed/DVD 2 

Flux (10“ 13 erg s _1 cm -2 ) 
AGN C model 3 

Component A Component B 

SBN F 6 

SBN E 7 

O VI A .1032+. 1037 

37.4/4334. /402. 

65.3+34.8 

6.6+3.5 


0.25+0.15 

(Lya) n A .1215 

101.8/3562./602 

701. 

1180. 



N IV] A .1487 

5.1/103./22.9 

21.4 

18.1 



(CIV)„ A .1549 

39.7/790. /177. 

384. 

718. 

... 


Hell A .1640 

21.4/426./95.5 

147. 

308. 

75. 

108. 

[Ne V] A .3426 

15.7/95./38.7 

115. 

152. 


10.9 

[Ne III] A .3869+. 3968 

19.2/97./43.2 

2.0+0. 6 

14.8+4.4 

42. +12. 

43. +13. 

Hell A .4686 

6.1/27.6/13. 

20. 

43.4 

11.3 

16. 

[O III] A .4959+.5007 

256./964./496 

13.4+4.6 

112. +39 

29. +83. 

156+452 

[Si VI] A 1.96 

8.0/9.2/8.6 

10.3 

8.7 



[Si VII] A 2.48 

8.3 

6.2 

0.9 



[Si IX] A 2.584 

3.0 

0.4 



... 

[Mg VIII] A 3.028 

11. 

3.6 




[Si IX] A 3.936 

5.4 

0.7 




[Mg IV] A 4.487 

7.6 

3.0 

13.3 


4.9 

[Ar VI] A 4.529 

15. 

18.9 

52.6 



[Mg VII] A 5.503 

13. 

11.2 

1.7 



[Mg V] A 5.610 

18. 

13.3 

34.5 


5.3 

[Ar II] A 6.985 

13. 



6.0 

2.0 

[Na III] A 7.318 

5.8 


... 

0.8 

0.5 

[Ne VI] A 7.652 

110. 

17L 

102. 


5.6 

[Fe VII] A 7.815 

3.0 

3.4 

9.2 



[Ar V] A 7.902 

< 12. 

1.7 

8.0 


1.2 

[Na VI] A 8.611 

< 16. 

1.4 

1.4 


0.25 

[Ar III] + [Mg VII] A 8.991 

25. 

••• +13.6 

0.8+2.2 

20. + 

7.3+ 

[Fe VII] A 9.527 

4.0 

13.7 

38.7 



[S IV] A 10.510 

58. 

13.5 

104 

1.3 

41.6 

[Ne II] A 12.813 

70. 



7.1 

2.5 

[Ar V] A 13.102 

< 16. 

1.8 

11.8 


2.2 

[Ne V] A 14.322 

97. 

133. 

355. 

... 

43.8 

[Ne III] A 15.555 

160. 

0.5 

5.6 

59. 

38. 

[S III] A 18.713 

40. 

0.5 

5.8 

21.3 

22.6 

[Ne V] A 24.317 

70. 

59.4 

280 


48. 

[O IV] A 25.890 

190. 

22. 

217. 

3.3 

207. 

[S III] A 33.481 

55. 

0.1 

3.3 

32.7 

34.8 

[Si II] A 34.814 

91. 

... 

0.1 

155. 

60. 

[Ne III] A 36.013 

18. 


0.5 

5.2 

3.4 

[O III] A 51.814 

no. 

0.2 

5.9 

16.7 

59. 

[N III] A 57.317 

51. 

0.1 

2.7 

5.9 

21.2 

[O I] A 63.184 

156. 



25.1 

8. 

[O III] A 88.356 

no. 


1.5 

26.6 

104. 

[N II] A 121.897 

30. 



7.0 

2.5 

[O I] A 145.525 

12. 



2.4 

0.8 

[C II] A 157.741 

220. 


... 

67.5 

38.5 


1 Dereddened line flux, assuming E b-v — 0.4 
2 Dereddened line flux, assuming Eb_v = 0-2 

3 AGN C parameters: same as AGN A models, but with the ionizing spectrum that includes a big blue bump 
(see text) 

4 SBN F parameters: Log U=-3.5, Log n=1.3, ionizing spectrum from Starburst99 continuous model with 
Z=.008 and duration of 50 Myr. The integration was stopped at a temperature of 50K. 

5 SBN E parameters; Log U=-2.5, Log n=0.3, ionizing spectrum from Starburst99 continuous model with 
Z=.008 and duration of 50 Myr. The integration was stopped at a temperature of 50K. 
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Table 4. Comparison of observed line fluxes with Ring Starburst model predictions 


Line id. A 
(/im ) 

Observed 

SBR A 1 

Flux (10 
SBR B 2 

13 erg s -1 
SBR C 3 

cm -2 ) 
SBR D 4 

SBR E 5 

SBR F 6 

[Ne II] A 12.813 

70. 

1.2 

1.5 

1.0 

1.2 

0.3 

1.1 

[Ar V] A 13.102 

< 16. 





0.3 


[Ne V] A 14.322 

97. 





0.6 


[Ne III] A 15.555 

160. 

11.2 

11.4 

8.7 

8.8 

5.7 

9.0 

[S III] A 18.713 

40. 

3.9 

3.7 

3.2 

3.0 

3.4 

3.2 

[Ne V] A 24.317 

70. 





7.3 


[O IV] A 25.890 

190. 



0.5 

0.4 

31. 

0.5 

[S III] A 33.481 

55. 

13.3 

12.9 

11.1 

11. 

12. 

11.3 

[Si II] A 34.814 

91. 

65.8 

75.2 

52.4 

58.4 

20.9 

53.8 

[Ne III] A 36.013 

18. 

2.2 

2.3 

1.8 

1.8 

1.1 

1.8 

[O III] A 51.814 

110. 

4.2 

3.7 

16.5 

14.7 

54.8 

15.3 

[N III] A 57.317 

51. 

1.5 

1.4 

5.9 

5.6 

19.5 

5.4 

[O I] A 63.184 

156. 

32. 

36.6 

21.7 

24.3 

7.3 

22. 

(O III] A 88.356 

110. 

7.0 

6.5 

26.6 

25.9 

95.8 

24.7 

[N II] A 121.897 

30. 

8.2 

8.0 

6.1 

6.2 

2.4 

6.7 

[O I] A 145.525 

12. 

3.1 

3.6 

2.1 

2.3 

0.7 

2.2 

[C II] A 157.741 

220. 

76.3 

142. 

60. 

107. 

35. 

62.3 


l SBR A parameters: Log U=-3.5, Log n=1.3, ionizing spectrum from Starburst99 instantaneous 
model with Z=.008 and duration of 5 Myr. The integration was stopped at a temperature of 50K. 

2 SBR B parameters: Log U=-3.5, Log n=0.3, ionizing spectrum from Star burst 99 instantaneous 
model with Z=.008 and duration of 5 Myr. The integration was stopped at a temperature of 50K. 

3 SBR C parameters: Log U=-3.5, Log n=1.3, ionizing spectrum from Starburst99 continuous 
model with Z=.008 and duration of 5 Myr. The integration was stopped at a temperature of 50K. 

4 SBR D parameters: Log U=-3.5, Log n=0.3, ionizing spectrum from Starburst99 continuous 
model with Z=.008 and duration of 5 Myr. The integration was stopped at a temperature of 50K. 

5 SBR E parameters: Log U=-2.5, Log n=0.3, ionizing spectrum from Starburst99 continuous 
model with Z=.008 and duration of 50 Myr. The integration was stopped at a temperature of 50K. 

6 SBR F parameters: Log U=-3.5, Log n=1.3, ionizing spectrum from Starburst99 continuous 
model with Z=.008 and duration of 50 Myr. The integration was stopped at a temperature of 50K. 
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Table 5. Comparison of observed line fluxes with composite model predictions 


Line id.A 
(*im ) 

Flux (10- 
Observed/DVD 2 

13 erg s~ 
CM1 3 

1 cm -2 ) 
CM2 4 

CM3 5 

O VI A .1032+. 1037 

37.4/4334. /402. 

170. 

26.4 

110.2 

(Lyo)n A .1215 

101. 8/3562. /602 

2426. 

2060. 

1881. 

N IV] A .1487 

5.1/103./22.9 

34. 

33.3 

39.5 

(CIV)„ A .1549 

39.7/790./177. 

774. 

512. 

1102. 

Hell A .1640 

21.4/426./95.5 

336. 

363. 

530. 

[Ne V] A .3426 

15.7/95./38.7 

137. 

70.1 

267. 

[Ne III] A .3869+. 3968 

19.2/97./43.2 

137.4 

180.3 

75.8 

Hell A .4686 

6.1/27.6/13. 

48. 

52.4 

74.7 

[O III] A .4959+. 5007 

256-/964. /496 

835 

889. 

281. . 

[Si VI] A 1.96 

8.0/9.2/8.6 

13.3 

14. 

19. 

[Si VII] A 2.48 

8.3 

5.5 

7.6 

7.1 

[Si IX] A 2.584 

3.0 

0.7 

1.0 

0.4 

[Mg VIII] A 3.028 

11. 

3.0 

4.6 

3.6 

[Si IX] A 3.936 

5.4 

1.2 

1.8 

0.7 

[Mg IV] A 4.487 

7.6 

15.6 

16.7 

16.3 

[Ar VI] A 4.529 

15. 

30.7 

31.4 

71.5 

[Mg VII] A 5.503 

13. 

8.7 

9.1 

12.9 

[Mg V] A 5.610 

18. 

25.4 

15.4 

47.8 

[Ar II] A 6.985 

13. 

7.4 

6.0 

6.0 

[Na III] A 7.318 

5.8 

1.4 

1.8 

0.8 

[Ne VI] A 7.652 

110. 

128.5 

71.1 

273. 

[Fe VII] A 7.815 

3.0 

7.5 

5.3 

12.6 

[Ar V] A 7.902 

< 12. 

5.3 

7.2 

9.7 

[Na VI] A 8.611 

< 16. 

1.5 

0.9 

2.8 

[Ar III] + [Mg VII] A 8.991 

25. 

42.6 

41.5 

36.6 

[Fe VII] A 9.527 

4.0 

31.4 

22.2 

52.4 

[S IV] A 10.510 

58. 

113. 

175.3 

118.8 

[Ne II] A 12.813 

70. 

14.7 

8.3 

7.4 

[Ar V] A 13.102 

< 16. 

7.4 

9.7 

13.9 

[Ne V] A 14.322 

97. 

230.6 

148.6 

488.6 

[Ne III] A 15.555 

160. 

115.1 

140. 

70.8 

[S III] A 18.713 

40. 

69.2 

55.7 

30.7 

[Ne V] A 24.317 

70. 

142.3 

85. 

346.7 

[O IV] A 25.890 

190. 

177. 

197.5 

273.3 

[S III] A 33.481 

55. 

67.7 

61. 

48.1 

[Si II] A 34.814 

91. 

177.4 

177. 

176. 

[Ne III] A 36.013 

18. 

10.5 

12.6 

6.8 

[O III] A 51.814 

110. 

115.5 

121.2 

77.6 

[N III] A 57.317 

51. 

40.6 

39. 

28.2 

[O I] A 63.184 

156. 

32.4 

32.4 

32.4 

[O III] A 88.356 

110. 

132.8 

134. 

123.9 

[N II] A 121.897 

30. 

9.6 

9.4 

9.4 

[O I] A 145.525 

12. 

3.1 

3.1 

3.1 

[C II] A 157.741 

220. 

103. 

102.5 

102.5 


J Dereddened line flux, assuming E b-v — 0.4 
2 Dereddened line flux, assuming E b-v — 0.2 
3 CM1 = AGN A + SBN F + SBR E 
4 CM2 = AGN B + SBN F + SBR E 
5 CM3 = AGN C + SBN F + SBR E 
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Table 6. Comparison of observed line fluxes with OH montecarlo model predictions 


Line id. A 


Flux (10 12 erg s 1 cm 2 ) 


Observed 


Modeled 
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